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The range over which thermodynamic temperatures have been realized by gas thermom- 
etry at the NBS has been extended to 730 K. The results are preserved by measuring the 
corresponding international practical temperatures. The difference between them is expressed 
as the following polynomial : 

7 1 /K-7 7 68 /K 68 =- 120,887.784/ T 6 i-f 1213.53295/ T 68 -4.3159552 + 6.44075647 X 10- 3 T 68 

- 3.56638846 X10" 6 T 6 2 8 

which is valid in the range 273 to 730 K. 

The difference found and the estimated uncertainties at the three defining fixed points 
in the range covered are 
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1 . Introduction 

The National Bureau of Standards' gas thermom- 
eter has been described in numerous papers [1-10] 1 
that give details of many of the important parts of 
the equipment and of some of the special measure- 
ments required. The last two of these papers also 
give the results of measurements with the gas ther- 
mometer from to 142°C. A comparison of these 
results with most earlier work elsewhere shows there 
are systematic differences that we ascribe to our more 
complete elimination of the effects of sorption. It was 
stated in the last paper (and has been further con- 
firmed by more recent observations) that our gas 
thermometer system and the helium filling gas were 
clean enough that any residual effect of sorption on 
these measurements would probably be negligibly 
small. Ultimately, it will be rewarding to remeasure 
these temperatures, when the benefits of further 



1 Figures in brackets indicate the literature references at the end of this paper. 



experiences with the gas thermometer can be realized. 
At present, it is more productive to extend the 
measurements to higher temperatures, both for their 
own intrinsic interest, and as a further step leading 
to the gold point. This paper gives the results of the 
initial measurements of the difference between ther- 
modynamic temperatures and international practical 
temperatures, at temperatures in the range from 141 
to 457 °C, and combines them with the earlier data 
[9, 10], reappraised in the light of further information. 

2. Equipment 

The essential elements of the gas thermometer 
have been described in earlier papers [9, 10]. Briefly, 
it is a constant volume t} r pe comprising a 432.5 cm 3 
bulb of platinum-rhodium alloy in the form of a right 
circular cylinder connected by small bore 90 percent 
platinum 10 percent rhodium tubing to a valve. The 
gas thermometer can be connected through the valve 
with a diaphragm pressure transducer, so that there- 
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after equality of the gas thermometer pressure with 
the pressure in the manometer can be established. 
The gas thermometer has a counter-pressure arrange- 
ment, a narrow annulus about 0.25 mm wide formed 
around the bulb by a heavy Inconel 2 600 case sur- 
rounding the bulb. 

The gas thermometer has been improved by in- 
stalling more thermocouples to measure the temper- 
ature distribution along the small bore tube, and by 
modifying the gas handling and vacuum equipment. 3 
The new arrangement of the vacuum system, shown 
in figure 1, allowed the system to be subdivided into 
four parts (labeled I, II, III, and IV) plus the 
manometer, so that all the critical sections could be 
pumped by ion pumps. The subsystems were ar- 
ranged in such a way as to be able to attain the best 
vacuum in part II, which included the gas thermom- 
eter bulb in direct communication with the residual 
gas analyzer (RGA). 

Two discrepancies from the earlier papers [9, 10] 
should be noted. The composition of the bulb is not 
all 80 percent Pt-20 percent Rh as stated in [9], nor 
all 88 percent Pt-12 percent Rh as stated in [10]; 
it is some of each. The sheet supplied for the top 
was found to have been 80 percent Pt-20 percent Rh, 
but the sheet for the sides and bottom was 88 



2 The use of commercial product names is for convenience and exactness and 
does not imply endorsement by the National Bureau of Standards. 

3 A new oil diffusion pump has been installed. It is stated by the manufacturer 
to have a backstreaming rate of "almost zero— less than 5X10 -5 mg per cm 2 
per min." 



percent Pt-12 percent Rh. However, the bulb is 
the same Bulb III referred to previously; the effect 
of this difference will be considered in the Discussion. 
Also, the volume of 454.8 cm 3 given for Bulb III 
has been re-evaluated as 432.5 cm 3 . The effect of this 
difference on earlier results will be tabulated in 
table 4, section 8. 

The liquids in the stirred baths were changed from 
organic materials. The liquid for the bath at °C 
was changed to water, to which potassium chromate 
was added as a rust inhibitor, and the liquid for 
the bath at temperatures from 142 to 457 °C was 
changed to a eutectic mixture of molten potassium 
nitrate, sodium nitrate and lithium nitrate. Despite 
the lesser viscosity of the new fluids, no important 
change in the performance of the thermostat was 
observed. 

New equipment was installed for the electrical 
measurements. All resistances of the platinum 
resistance thermometers (PRT's) were measured 
with improved speed and precision on an ac bridge 
[11]. All thermocouple emfs were measured by a 
digital voltmeter. 

For thermal expansion measurements, a new 
interferometer furnace, designed for use up to the 
gold point, was assembled, and it was operated from 
— 25 to 550 °C so that the thermal expansion of 
samples of the bulb material could be measured. 
Details of this equipment and the results will be 
given in a separate paper. 



Manometer 




Figure 1. Schematic of vacuum system. 
Boundaries of Parts I, II, III, and IV, which can be isolated for pumping, are shown by broken lines. 
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The apparatus shown in figure 2 was constructed 
for the determination of therm omolecular pressure 
effects. It consisted of two tubes, an 0.8 mm i.d., 
1.6 mm o.d. stainless steel tube running inside a 
9.6 mm i.d. Inconel tube, connected across a 
diaphragm pressure transducer. The outer tube was 
chosen of such a size that it could be inserted into 
one of the PRT wells of the gas thermometer bulb 
case to thermostat it in the stirred molten-salt bath. 
The apparatus was attached to the gas handling 
system, so that it could be evacuated, and thereafter 
gas could be introduced. The overall pressure was 
read with sufficient accuracy on a simple U-tube 
mercury manometer. The valves in the system were 
arranged so that when valves 2 and 4 were closed 
and 1 was open, the null of the diaphragm could be 
read, and when 1 and 3 were closed and 2 and 4 were 
opened, the thermomolecular pressure difference 
could be read on the diaphragm gage. 

3. Preliminary Considerations 

A study of the early results in the extended 
temperature range showed the need for modifying 
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Figure 2. Apparatus for measuring thermomolecular pressure. 



the experimental procedures. This section presents a 
discussion of some of those problems and their 
solutions. 

During vacuum bakeout we observed the ubiquity 
and persistence of hydrogen that is a common 
experience, and for which there has been a common 
explanation: that the hydrogen was absorbed in the 
metal in great quantity and diffused out of it slowly. 
Perhaps it was a happy accident that we introduced 
some wet helium into the counter-pressure system 
at 700 °C, and to our astonishment observed a 
drastic rise of the hydrogen peak on the RGA scan 
of the gas thermometer contents within 20 s. Inas- 
much as the helium peak did not change, the hy- 
drogen must have diffused through the wall, made of 
88 percent Pt-12 percent Rh metal. The diffusion of 
hydrogen through platinum is covered in the 
literature, which is ample, and is adequately sum- 
marized in Dushman [12]. There is considerable 
variation in the results, but a typical set of values at 
750 °C is as follows : 

Hydrogen diffuses through pure Pt at the rate of 
0.13 Pa-liter per cm 2 per min per mm thickness of 
the metal, with a hydrogen pressure of 10 5 Pa on 
one side and a vacuum on the other (the "permea- 
tion rate"). Platinum immersed in hydrogen at a 
pressure of 10 5 Pa absorbs 0.34 Pa-liters of the gas 
per cm 2 per mm thickness. 

Both the permeation rate and the solubility for 
hydrogen vary as P Vl , so their relative relationship 
remains the same at lower concentrations. The in- 
escapable conclusion is that hydrogen initially dis- 
solved in the platinum will quickly be exhausted, 
and converse^ that persistent concentrations of 
hydrogen observed during vacuum bakeout must be 
maintained by diffusion through the platinum. The 
same range of values of diffusion and permeation 
can be expected for platinum-rhodium alloys. The 
concentration of free hydrogen in the atmosphere is 
very small; the only possible sources are water or 
other hydrogen bearing materials which decompose 
to release hydrogen. Therefore, we surmised that 
the problem might be resolved by control of the 
ambient atmosphere. The furnace was flushed with 
argon during the bakeout, and a slow flow main- 
tained thereafter. Previously, the partial pressure of 
hydrogen in the gas thermometer bulb had never 
been lower than 70 /xPa, but it declined to 20 /xP& 
at 700 °C when the protective atmosphere was 
provided. 

For the present temperature range up to 457 °C, 
we found it necessary to change our procedures in 
three other ways which will be described in the 
following paragraphs : 

Modification oj procedure to avoid creep of the gas 
thermometer bulb. 

Above 327 °C, the gas thermometer bulb is sub- 
ject to creep (non-elastic deformation under small 
stress). It was discovered from the early results of 
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measurements carried out in the same manner as 
for the 0-142 °C range. The manometer equipment 
by its nature permits the accurate realization of 
only one pressure a day, with one or more experi- 
mentally determined values of temperature of the 
gas thermometer necessary to balance that pressure. 
Overnight, between measurements, we have always 
pumped out those portions of the system immediately 
adjacent to the gas thermometer bulb (parts II and 
IV) in order to maintain the purity of the system. 
(The pumpout of section II is also desirable for one 
of the "integrity checks' ' described below.) The 
typical pressure in the gas thermometer bulb at 
457 °C was 10 5 Pa; it was expected that it would 
deform elastically when the counter pressure annulus 
(in part IV) was evacuated. However, it was found 
that gas thermometer temperatures determined by 
measurements made in descending order agreed with 
those determined by measurements made in ascend- 
ing order only up to 327 °C. Above 327 °C, the gas 
thermometer temperatures differed more as the 
temperature (and, to some extent, time) increased. 
We also observed some hysteresis of the gas ther- 
mometer bulb after the counter-pressure was re- 
stored. It was necessary, therefore, in all subsequent 
measurements to keep the pressure inside and out- 
side the bulb nearly equal until the bulb temperature 
was below 300 °C Because it was more convenient 
to equalize the pressure by reducing it in the counter- 
pressure system, all values reported in the paper 
were derived from runs starting at the highest 
temperature and proceeding by consecutive steps of 
about 33 °C to the lower limit. 

Modification of Temperature Measurement of the 
Dead Space. 

The effect of the volume in the connecting tube 
of the gas thermometer, the "dead space," expressed 
in kelvins, is about six times as large at the zinc 
point as at the steam point. The uncertainties of 
the dead space effect can best be evaluated by 
examining the combined equation for the dead- 
space correction, 8t ds . This is 



&t ds =- 



zvr 2 /Vt i _p 1 Vt,\ 

V ^ k \T k , P.TJ 



(1) 



where the summation is over the length of the con- 
necting tube in mm increments, T± is the fiducial or 
reference temperature, T 2 is the measuring tempera- 
ture, Pj and P 2 are the corresponding pressures, V 
is the volume of the bulb at T u V ki and T ki are the 
corresponding values of volume and temperature 
at the &th position for state i. The element of volume 
varies between the two states only because of ther- 
mal expansion, i.e., V ki =V k (l+p ki 8t ki ), where p kt is 
the coefficient of thermal expansion depending upon 
the position and the state, and 8t ki is the difference 
in temperature of the kih position in the iih state 
from the calibrating temperature of 23 °C. The 
determination of the volume is thought to be suf- 



ficiently accurate [7]; the measurements of the T u >, 
on the other hand required improvement. 

We note from the equation that for any region of 
the dead-space actually at the gas thermometer 
temperature in both states, the net contribution 
is zero. The accuracy is increased, then, by main- 
taining the level of the liquid in the bath as high as 
possible. Unfortunately, the molten salt, because of 
low viscosity and low surface tension, could not be 
confined as readily as the earlier bath fluid. Even at 
the highest temperature the bath had to be less 
completely filled or -the salt would overflow, and it 
contracted so much over the range of temperature 
that the level was appreciably reduced at the lowest 
temperatures. 

Previously the temperatures for the dead space 
calculations for the steam point were derived on the 
following assumptions : 

(a) That the lowest junction is at the temperature 
of the bath. 

(b) That the emfs of the Pt-10 percent Rh/Pt 
thermocouples installed for dead space meas- 
urement follow the same temperature de- 
pendence as was determined for annealed 
thermocouples of the same spools of wire. 

(c) That the portion of the dead space above the 
uppermost junction was at the temperature 
of that junction. 

(d) That the values of the temperature can be 
adequately interpolated from least squares 
fits of the observed temperatures to poly- 
nomials chosen according to the complexity 
of the temperature profile. 

The consequence of these assumptions is that one 
may measure all emfs as differences from the bottom 
one and then the temperatures can be derived from 
the calibrations b}^ calculating the equivalent emf 
of each couple with an ice junction. 

Probably the steam point dead space calculation 
was adequately accurate because its contribution in 
kelvins was comparative^ small (0.034 K), the 
height of the bath liquid was satisfactory, the 
gradients were not severe, and the actual ther- 
mocouple emfs were reasonably close to standard 
table values. However, in this present temperature 
range, all these factors are less ideally fulfilled. 
The value of the dead space correction becomes large, 
partly because the bath liquid level had to be lower; 
the gradients, of course, were substantial, and as 
the temperature became high, the deviation of the 
actual emfs versus temperature from the table 
values became important. It is certain moreover 
that the bottom thermocouple was not adequately 
immersed. All these factors required that a better 
measurement and calculation be employed. 

Instead of referring the temperatures to the bottom 
thermocouple, if the temperature of the reference 
junctions of the dead space thermocouples were 
established, the error in the dead space calculation 
could be reduced. This is true because the temperature 
is then more accurately known than before at the 
lower temperatures, where the relative effect of a 
given volume of the dead space is larger; the same 
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discrepancy in emf produces a smaller error in 
kelvins if it affects the higher temperatures where the 
thermoelectric power becomes twice as large as at 
room temperature; and the accuracy of the result 
then no longer relies upon the bottom couple being 
at bath temperature. The details of this change will 
be discussed in section 6. 

Modification in the handling of the platinum resistance 
thermometers 

We found the platinum resistance thermometers to 
be much more sensitive to shock in the upper range 
of the temperature investigated. The vibration of 
the stirred liquid baths caused a significant increase 
of the resistances of the PRT* above 200 °C (not 
accountable by the effects reported by Berry [13], 
so that the procedures were changed to minimize the 
time during which the thermometers were kept in 
the bath. The measurements of the PRT resistances 
were made at the end of the time they were in the 
bath, therefore we used the values of the triple 
points that were found after each gas thermometer 
run to calculate their resistance ratios. With such 
precautions, the triple point values did not increase 
greatly, being subject to a change of <40 M^ from 
work hardening during a measurement period, and 
being partially restored to their original values by 
annealing at 425 and 457 °C whenever those tem- 
peratures were measured. 

4. Procedure 

The gas thermometer was prepared for operation 
by vacuum bakeout at 750 °C. After an initial 
pumpdown to a system pressure of 0.01 Pa by the oil 
diffusion pump, it was connected with the RGA and 
isolated from all other parts of the system. It was then 
pumped only by the ion pump of the RGA, a pro- 
cedure that gives both the highest concentration in 
the RGA for the detection of sorbable species, and 
very clean pumping for the gas thermometer. The 
entire bakeout furnace was isolated and a protective 
atmosphere of argon was introduced into it in order 
to exclude ambient air and prevent the diffusion of 
water vapor to the hot zone. The pumping was 
continued until the partial pressure of all contam- 
inating gases was less than 0.1 mPa in the gas 
thermometer. 

The measurements for this paper were made over a 
period of 15 months during which time the gas 
thermometer was loaded 10 times and measured at 
123 different temperatures. Of the first 51 tem- 
peratures measured only 9 aie free of creep effects, 
and of the remainder, the last 43 have the best 
measurements of the temperatures for the dead space. 

The gas thermometer was loaded while being 
thermostated at the highest temperatures to be 
measured. The initial measurements were commenced 
on the following day. At the beginning of a measure- 
ment, parts I, II, and III of the vacuum system were 
isolated but each was being evacuated, part I by its 
diffusion pump, part II by both ion pumps and part 



III b}^ its ion pump. The procedure for the measure- 
ment of a gas thermometer temperature was as 
follows : 

(1) Before loading, or at other times before 
introducing helium as a pressure transmitting gas, 
the diaphragm section was isolated by valves 1 
and 2, and the rate of pressure buildup was deter- 
mined over a period of a half hour or more. 4 The rate 
varied between the limits of 2.5 to 5 mPa/min. 
With a total volume of about 100 mm 3 , the relative 
amount of gas accumulating in the diaphragm com- 
pared with the amount of gas in the gas thermometer 
(for a pressure of 10 5 Pa) varied between 6 x 10~ 12 
min -1 and 1.2 x 10~ n min -1 . The measurement is an 
"integrity check" to ascertain that the diaphragm 
is neither contaminated nor has a leak from the out- 
side, and further that there is no significant leakage 
across the seat of valve 1. The same rate was ob- 
served whether the gas thermometer was evacuated 
or filled, hence it indicated no significant leakage 
occurred across the valve seat, and probably re- 
presented mostly, or totally, degassing. 

(2) At the same time, the heater for the Ti-CuO 
purification trap was turned on. Within 1 hr, the 
trap reached its operating temperature (700 °C), 
as confirmed by measurement with a Pt-10 percent 
Rh/Pt thermocouple referred to ice. 

(3) Two integrity checks of the ac bridge and 
triple point cell were made. First, the ratio of the 
standard 100 12 resistor to a 100 12 PRT was read 
to assure proper thermostating of the standard. 
Then the ratio of the resistance of a standard PRT 
immersed in a triple point cell to the standard resist- 
ance was determined on the ac bridge. This thermom- 
eter was treated with the utmost care, so that the 
ratio was reproducible within less than 1 [part in 10 7 
over a period of a month, a fact that tended to con- 
firm that the system was functioning properly. Over 
longer periods, the observed small changes in the ratio 
could have occurred because of drift in the resistances 
of either the standard resistance or the thermometer. 

(4) The ratios of the resistances of our 3 measur- 
ing PRT* to the standard resistor were determined 
at the temperature of the triple point of water. 

(5) The A4acLeod gage, monitoring the vacuum 
in the upper cell of the manometer was read (typi- 
cally 0.15 mPa or less) (an integrity check). 

(6) The ion pump currents were recorded (an 
integrity check). The RGA ion pump, open to part 
II, could be expected to have a current of 0.5 /xa, 
equivalent to 0.2 juPa pressure. 

(7) The Dewar around the molecular sieve (#13) 
trap was filled with liquid nitrogen, and when the 
temperature of the Ti-CuO trap had reached 700 °C, 
part I of the vacuum system was isolated and helium 
was slowly passed through the traps into the lines. 

(8) After isolating the ion pumps, helium was 
admitted to sections 2 and 4, until the pressure as 
indicated on the U-tube manometer was higher than 
the pressure in the manometer lines. Invariably, 
the gas flow was controlled to proceed in the direc- 



4 The subsequent discussion assumes the gas thermometer was already loaded. 
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tion from I to II and IV so that no backflow of 
helium would occur from other parts of the gas 
thermometer. Finally the manometer gas line was 
opened and the system pressure adjusted. When it 
was nearly at the pressure needed, the motor of the 
valve allowing communication of mercury between 
the upper cell and the lower cells of the manometer 
was turned on. The length of time it was necessary 
to open the valve before flow commenced was re- 
corded (an integrity check). Final adjustments of the 
amount of gas in the system were made, the mercury 
valve was opened fully and the level of mercury in 
the cells was checked. The mercury level was re- 
markably stable over periods of weeks (within about 
2 jum), a fact that depends not only upon a stable 
temperature but also upon a stable contact angle 
with the wall. After any needed adjustment of 
mercury level, the gas pressure was controlled by 
"thermal injection", governed by feedback from 
the capacitance bridge to maintain the pressure at 
the manometer setting. 

(9) After reading the diaphragm zero, we closed 
the by-pass valve and opened the valve that con- 
nected the gas thermometer valve with the dia- 
phragm. The pressure in the gas thermometer was 
made equal to the pressure in the manometer by 
adjusting the temperature of the liquid thermostat. 
When the pressures were about equal, the gas 
thermometer valve was closed, the by-pass valve 
opened, and the null of the diaphragm reread. The 
valves were then reset to reconnect the diaphragm 
to the gas thermometer so that any small difference 
between the pressure in the gas thermometer and the 
pressure in the manometer system could be observed 
and recorded. Then the gas thermometer was again 
isolated and the null reread. 

(10) The resistances of the PRT's in the gas 
thermometer bulb case were read, as well as that of 
the manometer reference PRT, then the 5 thermo- 
couples of the manometer and the 16 thermocouples 
of the dead space. 

(11) Another set of null-pressure-pressure-null 
measurements was made with the diaphragm. 

(12) Data also included were the time, the values 
of the gage blocks used in the manometer, and the 
regulator bridge and control settings. 

(13) The resistances of the platinum resistance 
thermometers were measured at the temperature of 
the triple point of water. 



5. Equations and Calculations 

The calculations for the gas thermometer tempera- 
tures can be made "symmetric", in that the equations 
can be devised so that the same calculation is per- 
formed for each measurement including that of the 
reference state. The details are given in [9]; the 
essentials are repeated here. We calculate a quantity 
designated Z that is an "augmented pressure", 
expressed in cm of Hg at 20 °C, where the thermal 
expansion, the pressure head and the dead space 
effects are accounted for. It is 



Z^hid+Wi), 



(2) 



where i refers to the iih measurement, h 0i is the 
height of the gage blocks used in the manometer, and 

+ [(T i +B i P i /R)/V ][j: k V k (l+t3 ki 8t ki )/T ki )]. (3) 

In the order of the terms in eq (3), the significance 
of the symbols is as follows : 



Ti is the pressure head and is given by 
*i=Ti*Mgl k /(RT kt ) 



(3a) 



where M is the molecular weight of the gas, g is the 
acceleration due to gravity, R is the molar gas 
constant=8.3137X10 6 cm 3 Pa mol" 1 K~\ T ki is the 
thermodynamic temperature for the kth element of 
length and ith measurement, and l k is the increment 
of length. 

The volume thermal expansion coefficient of the 
bulb, P, is derived from 

l+jK^a+aoOSandoo'^y 

as determined from samples of material of the bulb 
of length l at °C. The term 7r t /3 f is negligible. 

The remainder of eq (3) is the deadspace term. 
T t is the thermodynamic temperature of the bath, 
B t is the second virial coefficient of the thermometric 
gas at T f , P t is the pressure of the iih measurement, 
R is the molar gas constant, V is the volume of the 
gas thermometer bulb at °C, V k is the volume of 
the kth element of the connecting tube at 23 °C, /3 ki 
is the thermal expansion coefficient of the volume of 
the tube at temperature t kU in the kth position during 
the ^th measurement, and 5t ki =t ki —23 °C. This same 
temperature t ki , is expressed as an absolute thermo- 
dynamic temperature for T ki in the denominator. 
The quantity Z x is also given by the expression 



where 



Zi — ZjiQ-f-b/ji 

_ n R(T i +B i Pj/R) y 
PogV 



(4) 
(4a) 



with n the number of moles of gas, p the density of 
mercury at 20 °C, and other quantities as defined 
before. dZ t is a term for thermomolecular pressure, 
and any variation of the pressure head in the ma- 
nometer line from the level of the gas thermometer 
in the colder thermostat. We have calculated an 
"approximate gas thermometer temperature", T' GT , 
from the equation 



r' GT /iT GT =T ( /tf(Z,/Z 4 ). 6 



(5) 



Then the gas thermometer temperature, T GT , was 



* Scale comparisons will be given as pure numbers. All temperature scale quanti- 
ties are therefore divided by their corresponding unit. 
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derived by correcting T f GT for the effects of thermo- 
molecular pressures and differing pressure heads, 
so that 



REGION 



LENGTH 



Tgt/Kqt = T ' GT /K f GT + 5 1 P /K= 



■•Ti/KiZjJZn). 

(5a) 



We can substitute the expression for Z j0 and Z i0 
from 4a into 5a, and then because p , g, and V are 
constants, and n is maintained constant, we get 

T m /K m =T i /K+ ([Bj-BAPJPj) (Tj/T t )]P } /R)/K 

(6) 

with additional, insignificant second and higher order 
terms. Thus the gas thermometer temperature is 
equal to the thermodynamic temperature plus a 
term for the effects of gas imperfection. Substantially 

(P i /P j )(T j /T i ) = (l+kt j )=Con$L, for t j constant, 

(7) 

so that T GT /K GT =(T j +aP J )/K. Thus by deter- 
mining T GT /K GT at a constant pressure ratio but 
differing pressure ranges, Tj/K may be found as the 
intercept of the straight line of T GT /K GT versus P h 
and the difference of the second virial coefficients 
can be found from the slope. 

The calculation is made by computer in a series 
of programs: 

(1) The international temperatures of the gas 
thermometer and the manometer reference station 
are calculated by a program in Fortran IV named 
BRIDGE. This program is given in appendix I, 
with a file of the calibration constants A and B of 
the PRT's (TCAL75). A calibration on the ac 
bridge at the zinc, tin and triple points was made in 
June 1974, and the last gas thermometer measure- 
ments included in this paper were made in September 
of the same year. (The exact time does not appear to 
be important. Starting before the measurements 
were commenced, and extending to a time one year 
afterward, 4 sets of the thermometer calibrations 
have been shown to give a reproducibility of values 
at the temperatures of the fixed points within a 
standard deviation of 0.45 mK.) 

(2) The effect of the pressure head is calculated by 
a program in Basic named PYE. This program is 
given in appendix 2, and the subprogram LSFITl is 
given in appendix 3. In this program, the the treat- 
ment of the temperature distribution is made as 
described in section 6. The international tempera- 
tures, / 68 , were determined from the emf's of the 
thermocouples for the locations shown in figure 3. 
A correction was applied to / 68 to give the thermo- 
dynamic temperature, / th , and this in turn was 
converted to a value of reciprocal absolute tempera- 
ture, r. The volume of the gas in the constant volume 
valve and the first 6 mm of tubing within the valve 
were at the temperature of the valve. It is expected 




Figure 3. Thermocouple installation for measurement of the 
temperatures of the pressure head and dead space. 

The unit of length =0.929 mm. 



that the good thermal contact provided by the 
vacuum seal between the tube and the top of the 
header, H, brought the temperature of the tube 
substantially into agreement with the temperature 
of the header at that point. Between the valve and 
the top of the header (33 mm), the values of r were 
interpolated linearly, and similarly between the 
top of the header and the position of the thermo- 
couple Al (19 mm). The length of the tube above the 
top of the header varied by an amount D2, which 
reflected the differential expansion between the con- 
necting tube of the gas thermometer and its Inconel 
suspension. (Direct measurements confirmed the 
calculated value of D2 within experimental error.) 
The thermocouples Al . . ., A12, and B\ on the 
other hand terminated close to the connecting tube 
but were independently suspended, being installed 
in a harness that was self-supporting and through 
which the tube could move. The expansion of the 
thermocouple harness was included in the thermal 
expansion term for the tube. 

In order to interpolate the reciprocal temperature 
between the thermocouple locations, a quadratic 
equation, r=B(l)+5(2) /+2?(3) l\ was fitted by 
least squares to the first four pairs of values of r 
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versus position, /, for Al . . .A4 (76 mm). A 
polynomial giving r as 

r=B(l)// 2 +J3(2)//+5(3)+£(4) 7+5(5) I 2 

was fitted by least squares to the 9 pairs of reciprocal 
values of temperature versus position for ^44. . .^412 
(203 mm). In the interval from ^412 to B\ (25 mm) 
r was linearly interpolated. If the difference in emf 
between B\ and the calculated value for the tem- 
perature of the bath (the level of which was near or 
not more than 20 mm below Bl) was sufficiently 
large (>0.1 mV), the difference of the corresponding 
reciprocal absolute temperatures was linearly inter- 
polated over the next 136 units of length (126 mm), 
and if <0.1 mV, over 82 units of length (76 mm). The 
balance of the length of the connecting tube, and an 
additional length up to the center of the gas thermom- 
eter bulb, for a total length count of 688 (639 mm) 
was included in the calculation. The unit of length 
was 0.929 mm, which arose as a chart coordinate in 
measuring the tube diameter. The quantity w t is 
calculated as: 

as given in eq (3a), with the added definition that 
a ki is the linear expansion coefficient of the Mh 
element for the ?th measurement, and &t ki =t ki — 23. 
The output of the program gives the notebook 
reference number (S<t> is the notebook and page, Si 
the assigned run number), the value of 7r* in ppm, 
and the displacement D2 in mm. 

3. The calculation of the deadspace effect depends 
upon the evaluation of ^2 k Vki/T ki and involves 
exactly the same temperature calculations. They are 
associated with values of the volumes over the 
lengths from Jc-l to k. The calculation is performed by 
computer with a program in basic named DEDSPS, 
given in appendix 4. A file of chart readings for the 
diameter of the tube is a part of this program, that 
also required the entry of D2 from the PYE calcula- 
tion. The thermal expansion of the tube and of the 
thermocouple locations is combined in the term 
(l-\-fi ki At ki ). The output of the program gives the 
notebook and run reference numbers and the value 
of J^Vh/Tu, designated F*TAU. 

4. The value of Z is calculated by a program in 
Fortran IV named SUMMA, which is given in 
appendix 5. The calculation of the pressure is 
derived from eq (27) of [5], and is expressed in cm of 
mercury corrected to 20 °C. The output gives the 
notebook and run reference numbers and the value 
of Z in cm. 



6. Determination of Parameters 

Numerous calibrations and special measurements 
are involved in the gas thermometer calculations. 
Some of them have been mentioned in earlier 
sections, but they will be summarized here. 



6.1. The Gas Thermometer Pressure 

The value of pressure used in the gas thermometer 
equation is the sum of the pressure calculated by 
eq (27) of [5] plus the difference measured at the 
diaphragm. The height of the mercury column was 
determined by the length of the calibrated gage 
blocks. The density of the mercury was adjusted for its 
variation in temperature from 20 °C. We used the 
values of the thermal expansion of mercury pub- 
lished by Beattie et al. [14]. There is an offsetting 
effect of expansion of the gage blocks; the value in 
SUMMA reflects the difference between the volume 
expansion of mercury and the linear expansion of 
chromium carbide. The small values of the copper- 
constantan difference thermocouple readings for the 
manometer cells and for the main mercury column 
were converted to temperature differences by use of a 
"standard value" of the thermoelectric power, 
40.5 /xV/K [15], and the reference temperature was 
measured by a calibrated standard platinum re- 
sistance thermometer. 

In order to calculate the pressure differences from 
the gage readings of the diaphragm, its sensitivity 
was determined from the pressure changes calculated 
to have resulted from imposed changes of gas 
thermometer temperature. 

6.2. The Thermal Expansion of the Bulb 

The term in the gas thermometer equation second 
in importance to the pressure ratio is the term for 
the thermal expansion of the bulb. The thermal 
expansion coefficients of both 80 percent Pt-20 
percent Rh and 88 percent Pt-12 percent Rh have 
been carefully measured from —25 to 550 °C. The 
constants entered into SUMMA are those for the 
thermal expansion of 88 percent Pt-12 percent Rh 
where the length ratio is expressed as a polynomial, 



L/L =l + at+C 2 t 2 +C 3 t 3 +C^+C 5 t 5 . 



(8) 



The estimate of the standard deviation of a predicted 
point varies from 0.08 ppm in mid range to 0.14 ppm 
at either extreme. The estimate of the residual stand- 
ard deviation of 0.14 ppm is consistent with our 
expectation of the average imprecision of a single 
measurement. 

6.3. The Deadspace 

The values used in the deadspace calculations were 
obtained from various sources. The volume of the 
bulb was calculated from its dimensions except for 
the last gas load. In the bakeout preceding the final 
set of measurements, the bulb was deformed by an 
external overpressure. When the bulb was removed 
from the apparatus, its changed volume was deter- 
mined from the contained weight of water to be 
4.2103X10 5 mm 3 at °C, and this value was used 
in the calculations for the final set. 

To determine the volume of the tube, its diameter 
is required as a function of position. This was deter- 
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mined as described in [7] with an estimated total 
uncertainty that is insignificant in the final results. 
Temperatures were derived from thermocouple 
readings for the calculation of t and Vt. It may be 
useful to remark in advance that considerable error 
in the thermocouple values will not significantly 
affect the accuracy of the final results. The final 
arrangement illustrated in figure 3 was applicable 
to the last 30 states, where thermocouples were 
added about the header, II. The couples 66, 67, 
and (711 comprise copper-constantan legs (Type T), 
and 09 has Pt-10 percent Rh/Pt legs (Type S). 
The couple (76, with its reference junction at °C, 
was used to measure the temperature of the side of 
the header, and inferentially of the reference junc- 
tions of Al . . . A12, and J?l, because it is expected 
that the temperature of the reference ring D is nearly 
the same as the side of the header. The difference 
couples, (77 and (71 1 were used to determine the incre- 
ment of temperature between the top and side of 
the header and between the end of the tube at the 
constant volume valve and the side of the header, 
respectively. (79 was a difference couple running 
between the bath and the side of the header. It was 
not needed for the last 30 states, but was a necessary 
link with earlier measurements to evaluate the earlier 
data. All copper-constantan readings were converted 
to temperatures by use of a standard table [15]. The 
temperatures determined by all Pt-10 percent Rh/Pt 
thermocouples were evaluated from calibrations of 
similar couples which could be represented by a 
quartic, 

E=B l t+B 2 t+B 3 t 3 +Bj\ (9) 

with a standard deviation of the values at each inter- 
val of 25 °C up to 500 °C of 0.47 /xV when ^ = 5.45846 
X10- 3 M V/°C, 5 2 =1.13497X10" 5 /xV/(°C) 2 , B,= 
-1.52447X10- 8 mV/(°C) 3 and £ 4 = 9.06033X 10" 12 
M V/(°C) 4 . When the emf of the Pt-10 percent Rh/Pt 
thermocouple Bl was measured at 400 °C (under 
special conditions where the depth of immersion was 
thought to be adequate), however, it was found to 
differ from the equation by —28 /xV; the cause of the 
difference appeared to be work-hardening of the 
platinum leg during installation that was not subse- 
quently relieved b}^ annealing. This difference was 
assumed to be typical of all the deadspace couples 
and was taken into account in the gas thermometer 
calculations by multiplying all values of ^41 . . . 
A12, and 51 by 1.00895. 

In an earlier arrangment of the thermocouples 
applicable to the preceding 42 states, the (79 junction 
was located on the top of the header, and there was 
no (76, (77, nor (711 couple. Furthermore, (79 was 
calibrated at 400 ° C and found to differ from the 
value calculated from eq (9) by —51 /xV. Conse- 
quently, the measured emf's for (79 were multiplied 
by a factor, 1.01644, to account for the difference. 

To carry out the calculation described in the 
preceding paragraphs it was necessary to predict 
values of Al(5), 07(5), and 011(5), for the earlier 
measurements. This could be done from the infor- 
mation available from later measurements. In one 



switch position, the thermocouple emfs could be 
read with respect to J5l as the reference junction 
(position 4) and were designated as -41(4), etc. In 
another switch position, the thermocouples were ref- 
erenced to their own junctions with copper (position 
5), and measurements with this reference were 
designated as -41(5), etc. Because the temperature 
of the reference ring is uniform, .41(5)— -41(4) = 
51(5). Position 5 of the switch was used to read the 
thermocouples for the last 30 states. As shown in 
figure 4, the values of -41(5) and 67(5) are well- 
behaved functions of the bath temperatures as ex- 
pressed by (79 (5) . As for (71 1 (5) , its values depend not 
only upon the temperature of the bath, but also 
upon the temperature of the room and of the plate 
directly below it which is cooled by tap water. For 
constant room and cooling water temperatures, (711 
decreases linearly as the bath increases in tempera- 
ture, over a range of about 40 uV. The cooling water 
temperature varies seasonally from less than 8 °C 
in the winter to over 20 °C in the summer; this 
explains the large differences in (76(5). The 40 juV 
range of 611(5) tends to be more negative as the 
cooling water temperature increases. The values of 
611(5) used in table 2 were derived on the basis of a 
linear variation from 40 or 41 /zV at °C to txV 
at 457 °C. A uniform change of 20 juV in the values of 
(711(5) causers an error in the calculated gas ther- 
mometer temperature varying from 0.2 ppm at 140 
°C, up to 0.5 ppm at 457 °C. Therefore high accuracy 
in predicting these (71 1 (5) values is not a requirement. 




1000 2000 

C9(5) uY. 



3000 



Figure 4. Values of A 1(5) and C7(5) as a function of bath 
temperature, given by C9(5). 
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To summarize: For all the earlier data for which 
only values of (79(5) referenced to the top of the 
header are available, we find the temperature distri- 
bution along the tube by the following steps: 

(1) Kead (77(5) and Al (5) from figure 4. Calculate 
(711(5) by E-40 (457-^0/457 /*V. 

(2) Find the emf of a standard couple referred to 
ice for the top of the header in /xV for Pt— 10 percent 
Rh/Pt from V* (t tath )-E^X 1.01644, where V* is 
calculated from the quartic equation and E^ is (79(5). 

(3) Deduct (77(5) X 0.125 M V to account for the 
difference in emf corresponding to the difference in 
temperature between the top of the header and the 
aluminum reference ring. The equivalent emf for 
(76(5) can be found if desired. 

(4) Add Al (5) X 1.00895 to give a synthesized 
value of emf referenced to ice. 

(5) Find t A \ from the derived emf. Similarly, the 
temperatures can be found for all other positions 
including B\. 



6.4. Determination of Thermomolecular Pressures 

Using the apparatus developed for this purpose, 
and described in section 2, we measured some ther- 
momolecular pressure differences for helium at tem- 
peratures of 140, 260, 374, and 457 °C at pressures 
varying from 2.5X10 3 to >1X10 5 Pa. The results 
were correlated by an equation 



^<(WHW1 



(10) 



based on Weber's treatment [16], where the product 
of the average pressure p=(p 1 +p 2 )/2 and the pres- 
sure difference Ap is expressed in relation to T 2 , the 
upper temperature; T h the lower temperature (296 
K); and T , the reference temperature (273.15 K). 
The reference temperature is involved only in that 
the exponent n is defined by the assumption that the 
temperature dependence of the viscosity of helium 
can be expressed b} T 



^ =(r/T ) n+1/2 . 



(ii) 



For the temperature range °C <£ <460 °C, n 
was evaluated as 0.145. A least squares solution gave 
<7i = 3738 Pa 2 , with a standard deviation of 1.8 Pa 2 . 
The average deviation of the measured values from 
the equation is 1.8 percent, and the maximum devia- 
tion is 3.4 percent. Because C x varies inversely as the 
square of the diameter of the tube, the measured 
value of (7i was modified to account for the smaller 
diameter of the gas thermometer tube by a factor 
(0.823/0.904) 2 , or (7^3098 Pa 2 . 



7. Tabulation of Results 

The measurements which will be presented in sub- 
sequent tables are subdivided according to experi- 
mental groupings labeled Set A to Set I. More 
information by group is given in table 1. 





Table 1. Desig?iatio?is 


of sets, with added comme?its 


Set 


Period 


Pressure 
range 


Other 
comments 


A 

B 
C 
D 


Nov. 15-Dec. 5, 

1973. 
Jan. 15-Feb.l, 1974__ 
Feb. 6-Feb. 19, 1974_ 
Feb. 20-Mar. 5, 1974. 


Intermediate. _ 

Intermediate. _ 
Intermediate. _ 
High 


Effects of 
creep. 




Gas thermometer vacuum pumped at 457 °C May 10- 
May 15, 1974 


E 


May 16- June 6, 1974. 


High 






Resistance thermometers calibrated June 1-June 14, 
1974 


F 

G 


June 24- July 1, 1974 
July 6- Aug. 8, 1974__ 


Low 

Intermediate. _ 


Suspected con- 
tamination. 




Gas thermometer vacuum pumped at 457 °C 
August 9-August 12, 1974 


H 


Aug. 13- Aug. 30, 
1974. 


Intermediate. _ 






Gas thermometer vacuum Dumped at 920 °C 
August 31-September 11, 1974 


I 


Sept. 12-Sept. 13, 
1974. 


Intermediate. _ 


Change of bulb 
volume. 



The values of the calculated quantities for the 
states of the gas thermometer are tabulated in 
tables 2 and 3. The results have been divided into 
those in table 2, in which temperatures of the header 
were derived by calculation, and those in table 3 for 
which the thermocouple emf s were experimentally 
observed. The deviations of international tempera- 
tures from gas thermometer temperatures were cal- 
culated from the equation 

T ar/K' gt - T W /K M = TJKiZ/Zt) - T 68 /K 6S (12) 

where T t is the absolute temperature of the chosen 
reference state with the corresponding value of Z=Z t . 
The difference V GT !K' ' GT -T 68 /K &8 was further 
modified for the effect of the lower level of the gas 
thermometer bulb in the high temperature bath; the 
smaller pressure head requires a decrease in all 
the gas thermometer temperatures of 0.5 ppm except 



712 



Table 2. Gas thermometer results, sets A, B, C and D 



State 



2. 04 
11. 01 

1. 12 

1. 13 

6. 03 

8. 04 
10. 02 
10. 03 
10. 04 
10. 05 
10. 06 

9. 03 
9. 04 

8. 05 

10. 07 

9. 05 

8. 06 

7. 03 
6. 04 

6. 05 

5. 06 
4. 05 

3. 04 
1. 14 

1. 15 
3. 05 

2. 05 

11. 02 
1. 16 

10. 08 
10. 09 
10. 10 
10. 11 

9. 06 

8. 07 

7. 04 

6. 06 
1. 17 

12. 01 

12. 02 

13. 01 

14. 01 

15. 01 
1. 18 



Ave bath t 

°C 



199. 57969 
167. 81883 
. 05212 
. 00282 
327. 25501 
391. 69635 
457. 23666 
457. 23976 
457. 24714 
457. 23765 
457. 23818 
424. 93184 
424. 93024 
391. 71112 

457. 23859 
424. 92780 
391. 71202 
359. 55953 
327. 27613 
327. 27642 
295. 24589 
263. 30187 
231. 41455 
2. 40786 
— . 00855 
231. 41455 
199. 49531 
167. 82802 
. 08860 

457. 27365 
457. 27530 
457. 27417 
457. 27320 
424. 96377 
391. 74159 
359. 57478 
327. 29613 
. 02188 

400. 06997 
400. 07089 
323. 87915 
248. 07564 
172. 63305 
. 12744 
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8035 
1936 
5216 
271 
5619 
9901 
4089 
4431 
4748 
3943 
4196 
3321 
3153 
1223 

4473 

2889 

1381 

6005 

7628 

7648 

4676 

295 

1508 

787 

851 

1570 

9505 

2895 

8860 

7427 
7767 
7428 
7314 
6803 
4258 
7510 
9601 
2221 

7221 
7205 
7933 

7583 
3289 
2718 



703 



7965 
1882 
5187 
307 
5644 
9732 
3804 
4036 
5327 
3835 
4074 
3360 
3209 
1222 

3692 
2977 
1263 
6244 

7866 

7867 

4682 

171 

1502 

762 

874 

1449 

9614 

2784 

8849 

7405 
7593 
7685 
7602 
6335 
4223 
7576 
9679 
2182 

6941 
7250 
8026 
7615 
3369 
2756 



707 



7907 
1832 
5234 

268 
5239 
9272 
3107 
3459 
4066 
3518 
3185 
2870 
2709 

891 

3411 

2473 

961 

5610 

7344 

7409 

4409 

94 

1355 

810 

841 

1345 

9472 

2728 

8870 

7262 
7229 
7137 
7043 
5993 
3998 
7349 
9559 
2160 

6827 
6812 
7785 
7496 
3256 
2757 



Vault 

°C 



19. 9179 
19. 9868 
19. 9569 

19. 9797 

20. 0016 
20. 0017 
20. 0045 
20. 0285 
20. 0225 
20. 0272 
20. 0058 
20. 0049 
20. 0042 
20. 0034 

19. 9553 
19. 9956 
19. 9858 
19. 9158 
19. 9305 
19. 9298 
19. 9887 
19. 9734 
19. 9860 

19. 9939 

20. 0051 
20. 0092 
19. 9819 
19. 9700 
19. 9863 

19. 9872 
19. 9846 
19. 9864 
19. 9864 
19. 9716 

19. 9660 

20. 0132 

19. 9712 

20. 1028 

20. 0582 
20. 0345 
20 0446 
20. 0355 
20. 0243 
20. 0106 



C6 

mV 



575 
607 
611 
611 
547 
559 
587 
575 
560 
560 
500 
500 
500 
480 

667 
639 
635 
635 
627 
627 
667 
615 
579 
500 
524 
599 
587 
587 
555 

647 
647 
639 
639 
627 
587 
579 
575 
500 

619 
635 
615 
619 
595 
575 



CI 
mV 



19 
15 


35 
45 
54 
54 
54 
54 
54 
48 
48 
45 

51 
49 
45 
39 
36 
36 
31 
27 
23 


22 
18 
15 


54 
54 
54 
54 
49 
45 
40 
36 


45 
45 
36 
27 
17 
-1 



Cll 

mV 



26 

29 

41 

41 

14 

7 











3 

3 

7 


3 
7 
10 
14 
14 
17 
20 
23 
40 
40 
23 
26 
29 
40 









3 

7 

10 

14 

40 

6 
6 
14 
22 
29 
36 



10 



7. 60 

8. 03 
10. 95 
10. 96 

6. 51 
6. 00 
5. 67 
5. 67 
5.68 
5. 67 
5. 68 
5. 88 

5. 87 

6. 07 

5. 65 

5. 84 

6. 05 
6. 26 
6. 50 
6. 50 

6. 74 

7. 02 

7. 33 
10. 91 
10. 97 

7.35 
7.72 

8. 02 
10. 10 

5. 67 
5. 67 
5. 67 
5. 67 

5. 85 

6. 06 
6. 28 
6. 51 

10. 98 

6. 04 

6. 02 

6. 55 

7. 17 
7. 97 

10. 95 



11 



Vt 

mm?IK 



0. 843234 
. 879305 

1. 143320 
1. 144712 

. 730605 
. 677523 
. 644306 
. 644577 
. 645384 
. 644916 
. 646497 
. 667050 
. 665542 
. 687159 

. 642461 
. 661441 
. 683365 
. 704933 
. 729316 
. 729117 
. 752268 
. 780777 
. 811778 
1. 14086 
1. 14563 
. 814564 
. 851899 
. 879258 
1. 14442 

. 644047 
. 644099 
. 644346 
. 644343 
. 663317 
. 685320 
. 706876 
. 730882 
1. 14667 

. 683636 
. 681597 
. 734768 
. 796473 
. 874040 
1. 14332 



12 



Gage block 
in 



19. 6000184 

18. 3000156 
11. 3926111 
11. 3906113 
24. 8000249 

27. 4080254 
30. 0483289 
30. 0483289 
30. 0483289 
30. 0483289 
30. 0483289 

28. 7480296 
28. 7480296 

27. 4080254 

30. 0483289 

28. 7480296 

27. 4080254 

26. 1080274 
24. 8000249 
24. 8000249 

23. 5000240 
22. 2000212 

20. 9000191 
11. 4900096 
11. 3900086 
20. 9000191 

19. 6000184 
18. 3000156 
11. 3940123 

30. 0483289 
30. 0483289 
30. 0483289 
30. 0483289 

28. 7480296 

27. 4080254 
26. 1080274 

24. 8000249 
11. 3910128 

36. 0000341 
36. 0000341 
32. 0000308 

28. 0000326 
24. 0000336 
14. 7860234 



13 



Z 

cm 



50. 098428 
46. 732459 
28. 955979 

28. 950756 
63. 628995 
70. 456473 
77. 400202 
77. 400044 
77. 400303 
77. 399842 
77. 400268 
73. 977531 
73. 977283 
70. 457598 

77. 400091 
73. 976691 
70. 457391 
67. 050862 
63. 629693 
63. 629769 

60. 235626 
56. 850078 
53. 470746 

29. 205208 
28. 949104 
53. 470813 
50. 098911 
46. 723628 
28. 959419 

77. 400303 
77. 400310 
77. 400339 
77. 400290 
73. 977511 

70. 457874 
67. 049888 
63. 629483 
28. 951180 

92.567811 
92. 567881 
82. 093154 

71. 669921 

61. 296616 
37. 580545 



14 



t/Kc 



-0. 0468 
-. 0442 
-. 00002 
Fiducial 
-. 0601 
-. 0836 
-. 1092 
-. 1138 
-. 1188 
-. 1136 
-. 1101 
-. 0976 
-. 0984 
-. 0878 

-. 1009 
-. 0907 
-. 0803 
-. 0692 
-. 0654 
-. 0649 
-. 0589 
-. 0583 
-. 0558 
-. 00001 
Fiducial 
-. 0551 
-. 0506 
-. 0449 
. 00017 

-. 1050 
-. 1066 
-. 1052 
-. 1047 
-. 0913 
-. 0790 
-. 0686 
-. 0635 
Fiducial 

-.0873 

-. 0877 
-. 0658 
-. 0578 
-. 0479 
Fiducial 



15 



Tqt/Kgt — 



-0. 0462 
-. 0438 
-. 00002 



, 0587 

, 0821 

, 1075 

, 1121 

, 1170 

, 1119 

, 1084 

0960 

0967 

0863 

0992 
0891 
0788 
0679 
0641 
0637 
0578 
0573 
0549 
00001 



0542 
0500 
0446 
00017 

1033 
1049 
1034 
1029 
0896 
0775 
0672 
0622 



0865 
0869 
0653 
0571 
0473 



16 



TjK-T 6S /K ti 



-0. 0390 
-. 0381 
-. 0002 



-. 0444 
-. 0634 
-. 0842 
-. 0888 
-. 0937 
-. 0886 
-. 0851 
-. 0751 
-. 0758 
-. 0676 

-. 0759 
-. 0682 
-. 0601 
-. 0514 
-. 0498 
-. 0494 
-. 0454 
-. 0468 
-. 0461 
+. 00003 



0454 
0423 
0389 
00017 

0800 
0816 
0801 
0796 
0687 
0588 
0507 
0479 



-.0617 
-. 0621 
-. 0470 
-. 0444 
-. 0403 




Identification of columns: 1 
determined by 
sure head effect 
ence between 
momolecular pre 



i— l 













Table 3 


. Gas thermometer resui 


ts, sets E, F 


G, H and I 










1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


State 


Ave Bath t 


701 


703 


707 


Vault 


C6 


C7 


Cll 


7T 


Vr 


Gage blocks 


Z 


T'gt/Kqt — 


Tqt/Kgt — 


T/K- 




°C 








°C 


*v 


f*V 


nV 




mm*/K 


in 


cm 


T$s/Kqs 


Tss/Kss 


TwIKw 


9. 04 


400. 07781 


7938 


7913 


7493 


19. 9122 


1049 


50 


3 


5. 84 


0. 656597 


36. 0000180 


92. 566491 


-0. 0863 


-0. 0855 


-0. 0607 


9.05 


400. 07842 


8163 


7785 


7577 


19. 9145 


1041 


50 


14 


5. 84 


. 656825 


36. 0000180 


92. 566547 


-. 0882 


-. 0874 


-. 0626 


8. 01 


323. 88316 


8329 


8433 


8185 


19. 9505 


1007 


39 


27 


6. 37 


. 710510 


32. 0000120 


82. 091827 


-. 0620 


-. 0614 


-. 0431 


7. 01 


248. 08543 


8546 


8605 


8478 


19. 9266 


989 


27 


32 


7. 00 


. 773604 


28. 0000080 


71. 669463 


-. 0555 


-. 0551 


-. 0424 


6. 01 


172. 64554 


4527 


4590 


4546 


19. 8996 


921 


20 


45 


7.82 


. 853059 


24. 0000090 


61. 296727 


-. 0458 


-. 0455 


-. 0380 


1. 01 


. 21285 


1298 


1260 


1298 


19. 4946 


841 


4 


79 


10. 89 


1. 13583 


14. 7900210 


37. 591184 


-. 00009 


-. 00009 


-. 00009 


1. 02 


. 02653 
140. 68118 


2657 
8090 


2662 
8143 


2641 
8121 


19. 9315 
19. 9331 


853 
883 


5 
21 


82 
67 


10.90 
8. 23 


1. 13647 

. 893283 


14. 7800240 
22. 3000130 


37. 565563 
56. 902358 


Fiducial 
-. 0373 






5. 01 


"-."0371" 


""-."63I3" 


10. 01 


457. 29908 


9784 


xOOO 


9939 


19. 9697 


1153 


60 


2 


5. 49 


. 620079 


10. 5510989 


27. 177196 


-. 1042 


-. 0872 


-. 0790 


10. 02 


457. 29983 


9816 


x333 


9801 


19. 9738 


1158 


62 


-1 


5. 49 


. 619983 


10. 5510989 


27. 177176 


-. 1055 


-. 0885 


-. 0803 


11. 01 


379. 56253 


6289 


6353 


6117 


19. 9944 


1096 


49 


16 


5.97 


. 670790 


9. 4508962 


24. 285564 


-. 0757 


-. 0617 


-. 0555 


11. 02 


379. 56304 


6331 


6412 


6170 


19. 9939 


1093 


49 


20 


5. 96 


. 668803 


9. 4508962 


24. 285539 


-. 0769 


-. 0629 


-. 0567 


1. 03 


. 06066 
457. 29595 


6070 
9467 


6050 

9880 


6078 
9437 


19. 9752 

20. 0108 


914 
1318 


3 

57 


62 
-25 


10. 88 
5. 47 


1. 13410 
. 617593 


4. 0000033 
30. 0483217 


10. 166566 
77. 396673 


Fiducial 
-. 0960 






10. 03 


-. 0889 


-. 0656 


10. 04 


457. 29569 


9559 


9732 


9417 


20. 0127 


1318 


58 


-19 


5. 46 


. 616579 


30. 0483217 


77. 396573 


-. 0913 


-. 0896 


-. 0663 


10. 05 


457. 29379 


9446 


9548 


9143 


20. 0123 


1321 


55 


-46 


5. 47 


. 617047 


30. 0483217 


77. 396510 


— . 0900 


-. 0883 


-. 0660 


12. 01 


424. 97992 


8277 


8039 


7662 


20. 0264 


1302 


52 


-49 


5.66 


. 636436 


28. 7480235 


73. 973390 


-. 0785 


-. 0769 


-. 0560 


13. 01 


391. 76711 


6662 


6968 


6503 


20. 0098 


1231 


46 


-26 


5. 86 


. 657983 


27. 4080211 


70. 454443 


-. 0723 


-. 0708 


-. 0521 


14. 01 


359. 61197 


1100 


1368 


1122 


19. 9801 


1202 


40 


-5 


6. 08 


. 679571 


26. 1080190 


67. 047804 


-. 0640 


-. 0626 


-. 0460 


15. 01 


327. 32650 


2646 


2677 


2627 


19. 9761 


1228 


39 


-16 


6. 30 


. 701470 


24. 8000160 


63. 626934 


-. 0596 


-. 0584 


-. 0441 


16. 01 


231. 46721 


6677 


6738 


6747 


19. 9887 


1128 


24 


-13 


7. 15 


. 786948 


20. 9000176 


53. 469426 


-. 0519 


-. 0510 


-. 0422 


16. 02 


231. 46773 


6751 


6776 


6792 


19. 9878 


1129 


24 


-1 


7. 15 


. 786771 


20. 9000176 


53. 469488 


-. 0518 


-. 0509 


-. 0421 


17. 01 


136. 16997 


6928 


7004 


7058 


19. 9771 


1085 


12 


12 


8. 26 


. 893907 


17. 0000212 


43. 372362 


-. 0358 


-. 0355 


-. 0311 


1. 04 


. 02701 
457. 29574 


2679 
9756 


2715 
9418 


2709 
9549 


19. 9949 
19. 9837 


1040 
1165 


-1 
51 


24 
-11 


10. 85 
5. 51 


1. 12928 
. 623306 


11. 3900096 
30. 0483217 


28. 948912 
77. 397914 


Fiducial 
-. 1017 






10. 06 


-. 1000 


-. 0767 


10. 07 


457. 29670 


9884 


9768 


9358 


19. 9867 


1176 


52 


-12 


5. 51 


. 622538 


30. 0483217 


77. 397748 


-. 1043 


-. 1025 


-. 0792 


10. 08 


457. 29657 


9844 


9777 


9350 


19. 9857 


1179 


52 


-11 


5. 50 


. 622227 


30. 0483217 


77. 397772 


-. 1039 


-. 1022 


-. 0789 


12. 02 


424. 97991 


8042 


8006 


7926 


19. 9913 


1213 


46 


-27 


5. 68 


. 639977 


28. 7480235 


73. 974725 


-. 0879 


-. 0863 


-. 0654 


13. 02 


391. 76040 


5917 


6261 


5943 


19. 9937 


1190 


41 


-28 


5. 88 


. 660398 


27. 4080291 


70. 455153 


-. 0798 


-. 0783 


-. 0596 


14. 02 


359. 60084 


1 


126 


124 


19. 9980 


1171 


37 


-22 


6. 09 


. 682020 


26. 1080270 


67. 048033 


-. 0706 


-. 0693 


-. 0528 


15. 02 


327. 31830 


1827 


1782 


1882 


20. 0035 


1172 


30 


-20 


6. 36 


. 709625 


24. 8000240 


63. 627550 


-. 0645 


-. 0632 


-. 0490 


18. 01 


295. 29733 


9749 


9761 


9689 


19. 9961 


1154 


27 


-21 


6. 61 


. 734450 


23. 5000305 


60. 234739 


-. 0589 


-. 0577 


-. 0453 


19. 01 


263. 34529 


4570 


4505 


4513 


19. 9825 


1175 


22 


-20 


6.87 


. 760698 


22. 2000277 


56. 848808 


-. 0572 


-. 0562 


-. 0456 


20. 01 


231. 45545 


5538 


5528 


5568 


19. 9956 


1155 


20 


-6 


7. 17 


. 789847 


20. 9000256 


53. 469639 


-. 0540 


-. 0531 


-. 0443 


21. 01 


199. 62867 


2848 


2866 


2886 


20. 0017 


1156 


16 


-7 


7. 50 


. 821595 


19. 6000301 


50. 097231 


-. 0500 


-. 0494 


-. 0417 


1. 05 


. 01667 
199. 62796 


1671 
2762 


1654 
2790 


1677 
2836 


20. 0060 
19. 9792 


1105 
1075 


-3 
23 


3 

-9 


10. 83 

7. 48 


1. 12683 
. 819122 


11. 3900096 
19. 6000301 


28. 948732 
50. 098373 


Fiducial 
-. 0503 






21. 02 


-. 0497 


-. 0420 


1.06 


. 01724 


1734 


1727 


1711 


19. 9817 


990 


-3 


37 


10. 86 


1. 13133 


11. 3900096 


28. 949522 


Fiducial 























T\ 


\BLE 4. 


Gas thermometer results previously published [10] 










1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


State 


Ave bath t 

°C 


701 


703 


707 


Vault 
°C 


IT 


Vr 

mm 3 /K 


Gage blocks 
in 


Z 

cm 


T" gt/Kgt — 


5£buib K 


S^therm -& 


8t xmv K 


T/K- 


2. 02 
1. 03 


100. 00097 
-. 00468 

99. 99979 
-. 00460 

-. 00066 
100. 00427 
120. 78497 
141. 59376 
100. 00536 

10. 000484 

120. 79043 

141. 59935 

100. 00710 

. 00254 

20. 65967 

41. 35574 

62. 06274 

82. 79370 

100. 01211 

. 00580 

20. 66189 
41. 36032 
62. 06529 
62. 06526 

82. 78611 
. 00497 


0093 
464 

9972 
464 

59 

423 

8478 

9364 

515 

473 

935 
9925 
698 
252 
5969 
5582 
6282 

9344 

1195 

574 

6184 
6030 
6521 
6537 
8611 
500 


0065 
485 

9995 

477 

86 

397 

8464 

9334 

511 

447 
9011 
9898 
681 
243 
5955 
5556 
6257 

9333 

1186 

571 

6190 
6016 
6519 
6521 
8600 
479 


0133 
455 

x002 
441 

55 
461 

8548 

9429 

581 

531 
9084 
9983 
749 
269 
5976 
5583 
6283 

9434 

1253 

595 

6193 
6051 
6546 
6520 
8621 
514 


20. 0101 
20. 0128 

20. 0381 
20. 0239 

20. 1096 
20. 0166 
20. 0204 
20. 0146 

19. 9973 

20. 0091 
20. 0052 
20. 0029 
20. 0080 

19. 9962 

20. 0006 
19. 9952 
19. 9905 

19. 9716 
19. 9975 
19. 9929 

19. 9816 

19. 9789 

20. 0144 
20. 0145 
20. 0119 
19. 9822 


8.63 
10. 64 

8.62 
10. 64 

10. 64 

8. 64 
8. 30 
8. 03 
8.60 

8. 59 
8. 30 
8. 03 

8. 59 
10. 61 
10. 10 

9. 64 
9. 21 

8. 90 

8. 59 
10. 61 

10. 10 

9. 67 
9. 24 
9. 24 
8. 86 

10. 61 


0. 723352 
. 836102 

. 721825 
. 835937 

. 835661 
. 724165 
. 701367 
. 685585 
. 719052 

. 718361 
. 701608 
. 685669 
. 717602 
. 831801 
. 804239 
. 779314 
. 753338 

. 738540 
. 718271 
. 832181 

. 804688 
. 783029 
. 756769 
. 756665 
. 733769 
. 832388 


30. 0496312 
22. 0592231 

30. 0496312 
22. 0592231 

3. 9999996 
5. 4888939 

5. 7488897 

6. 0488946 
5. 4488939 

5. 4489006 

5. 7488981 

6. 0488977 
5. 4489006 

4. 0000047 
4. 3000022 
4. 6002021 

4. 9001996 

5. 2000040 
5. 4889006 
4. 0000047 

4. 3000022 
4. 6002021 
4. 9001996 

4. 9001996 

5. 2000040 
4. 0000047 


76. 569444 
56. 053471 

76. 568881 
56. 053319 

10. 164256 

13. 884388 

14. 567342 

15. 431416 

13. 884378 

13. 884336 

14. 657444 

15. 431489 
13. 884351 
10. 164242 

10. 932764 

11. 702679 

12. 472959 

13. 243896 
13. 884349 
10. 164321 

10. 932783 

11. 702738 

12. 472967 

12. 472965 

13. 243744 
10. 164280 


-0. 0325 
Fiducial 

-. 0330 

Fiducial 

Fiducial 
-. 0319 
-.0406 
-. 0473 
-. 0333 

-. 0290 
-. 0382 
-. 0455 
-. 0305 
Fiducial 
-. 0040 
-. 0094 
. 0155 

-. 0275 
-. 0344 
Fiducial 

-. 0045 
-. 0111 
-. 0166 
-. 0167 
-. 0236 
Fiducial 


+ 0.0019 


0. 0001 


0. 0058 


-0. 0247 


2. 03 
1. 04 


. 0019 


. 0001 


. 0058 


-. 0252 


1. 05 










2. 04 

3. 01 

4. 01 
2. 05 

2. 06 

3. 02 

4. 02 
2. 07 
1. 06 


.0019 
.0023 
. 0028 
. 0019 

. 0019 
. 0023 
. 0028 
. 0019 


.0041 
. 0050 
. 0057 
. 0041 

. 0041 
. 0050 
. 0057 
. 0041 


.0010 
. 0013 
. 0016 
. 0010 

. 0010 
. 0013 
. 0016 
. 0010 


-. 0249 
-. 0320 
-. 0372 
-. 0263 

-. 0220 
-. 0296 
-. 0354 
-. 0235 


5. 01 

6. 01 

7. 01 

8. 01 
2. 08 
1. 07 


. 0003 
. 0007 
. 0010 

. 0015 
. 0019 


. 0010 
. 0018 
. 0026 

. 00035 
. 0041 


. 0002 
. 0004 
. 006 

. 0008 
. 0010 


-. 0025 
-. 0065 
-. 0113 

-. 0217 
-. 0274 


5. 02 

6. 02 

7. 02 

7. 03 

8. 03 
1. 08 


. 0003 
. 0007 
. 0010 
. 0010 
. 0015 


. 0010 
. 0018 
. 0026 
. 0026 
. 0035 


. 0002 
. 0004 
. 0006 
. 0006 
. 0008 


-. 0030 
-. 0082 
-. 0124 
-. 0125 

-. 0178 















Identification of columns: Columns 1-6 as for table 2, 7 = relative pressure head effect, 8 = value for dead space £*Vjk/T*,-, 9 = gage block height used in manom- 
eter, 10 = Z= augmented pressure expressed in cm of Hg, ll = difference of values found from gas thermometer and PRT's, 12 = correction required for the revised 
volume of the bulb, 13 = correction for the thermomolecular pressure effect, 14 = correction for gas imperfection, 15 = difference of thermodynamic and International 
Practical Kelvin temperatures at the temperature of Us, given in column 2. 



those measured in the low temperature bath. These 
changes are smaller than the thermomolecular 
pressure effects that were calculated by eq (10), 
and that were included in the results in column 15, 
headed T GT /K GT -T, 8 /K 6S . 

8. Review of Previously Published Results 

The results reported in our last paper [10] were 
calculated by programs similar to those used for the 
present data. The PYE and DEDSPS programs 
derived temperatures in the manner described as 
our "original method," and the SUMMA program 
differed in the form of the equation and the constants 
used for the thermal expansion of the bulb, with a 
trivial effect on the results. Both of these differences 
can be expected to produce a difference less than 0.5 
mK from the present approach. The original data 
and some of the quantities derived from intermediate 
calculations are given in table 4, columns (1) to (10), 
and the values reported previously are given in 
column (1 1) . As stated in section 2, the actual volume 
of the bulb was found to be different from the value 
used in the original SUMMA program. The change 
necessitated by correcting the volume of the bulb 
is given in column (12). A further adjustment for 
the effects of thermomolecular pressure, as deter- 
mined by eq (10), is given in column (13). 

9. The Effect of Gas Imperfection 

When the numerical differences, T GT jK GT — 
T 68 /K &8 , for the same temperature obtained at nearly 
constant pressure ratio but over different ranges of 
pressure are extrapolated vs the upper pressure 
to zero pressure, the intercept is TjK — T &8 /K 68 . 
Alternatively, if the equation of state for 1 mol of gas 
is considered in the form pV=RT-{-Bp, the thermo- 
dynamic temperature T is related to a gas ther- 
mometer temperature T GT by 

T/K= T GT /K GT - { [B, -Btfa/p,) (Tj/TMpjW/K 

(13) 

where j's refer to the measuring state, and i's to the 
reference state, and R to the molar gas constant. 
Considered in terms of this equation, the accuracy 
with which the thermodynamic temperature can be 
found by extrapolation will be optimized when the 
range of p/s is large, so long as the lowest p t is still 
large enough to be realized with the same relative 
accuracy as large values. However, at low pressures, 
thermomolecular pressure effects become important. 
For instance, two pressure ranges were used at 
100 °C. For p, = lX10 5 Pa, the combined thermo- 
molecular pressure effects for p } and p x are 0.34 ppm, 
while for ^=18374 Pa, the combined thermo- 
molecular pressure effects are 10.6 ppm. For all 
values in Sets A, B, C, G, H, and I, all p/s were 
above 5X10 4 Pa, so that the thermomolecular pres- 
sure affects the gas thermometer temperature by at 
most 2 mK. 



Unfortunately, the number of measurements at 
varying pressure ranges was insufficient to permit 
precise extrapolation to zero pressure. An alternative 
to extrapolating is to calculate the effects of gas 
imperfection from independently determined second 
virial coefficients of helium. Our original decision to 
rely upon extrapolation involved our appraisal that 
the use of virials from the literature would not give as 
accurate thermodynamic temperatures. This re- 
mains our conclusion, provided there are sufficient 
gas thermometer measurements ; 6 in the interim we 
are compelled to use literature values of virials to 
derive thermodynamic temperatures, but with more 
uncertainty than we would expect when our own 
measurements are sufficiently complete. 

Second virial coefficients calculated from com- 
pressibility measurements are probably more accur- 
ate than from other sources; the Burnett method 
[17] is probably the most accurate of the experimental 
compressibility techniques. We have chosen a set of 
virial coefficients based upon the values obtained at 
the NBS by the Burnett method for temperatures up 
to 150 °C [18]. The values at higher temperatures 
are those of Yntema and Schneider [19], who also 
used the Burnett method. The interpolation and 
smoothing were made consistent with acoustic 
results of Gammon and Douslin [20]. The results 
from these three sources differ from each other sub- 
stantially more than their respective estimated 
uncertainties. The values used and the equation 
derived by a least squares fit to the data are given in 
table 5. 

Table 5. The second virial coefficient of Helium 



t(°C) 


B (cm 3 /mol) 





12. 00 


25 


11.89 


50 


11. 77 


75 


11.67 


100 


11. 56 


125 


11.46 


150 


11. 36 


300 


10. 76 


400 


10. 45 


500 


10. 14 


600 


9. 82 



B/(cm 3 mol) = 11.9967—4.48574 

X10- 3 £ + 1.46724X10~ 6 £ 2 . (14) 

The effects of gas imperfection increase with in- 
creasing temperature and, of course, with increasing 
pressure. It amounts to 0.0227 K for 2*2=730 K and 
p = \0 5 Pa, and only 0.0002 K for T 2 = 293.81 K 
and #>= 14500 Pa. The multiplying factor on B u 
f=pJpj.Tj/Ti, is always greater than unity above 
°C, therefore the sign of the difference of the 



6 However, another approach using low pressures corrected by carefully deter- 
mined thermomolecular pressures is probably better. It is discussed in the final 
section of this paper. 
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virials B 2 —jB X) is always negative. Consequently, 
the thermodynamic temperature calculated from 
eq (13) is always greater than the gas thermometer 
temperature. 

The results corrected for gas imperfection are 
given in the final column of tables 2, 3, and 4, and 
represent the difference between international and 
thermodynamic temperatures. An equation giving 
the deviation as 

T/K-Tss/K^Ar/T^ 

+A 2 /T 6S +A 3 +AJ\ s +A 5 T m > (15) 

was fitted to all the values in tables 2, 3, and 4 by 
least squares. The constants found are given in 
table 6. 

Table 6. — Constants for calculation of T/K—Tm/Km 
(eq (15)) 



The estimate of the standard deviation of the fit, 



A x = - 1. 208877838 X 10 5 AV 
A 2 = + 1.2135329499 X 10 3 K 68 
A 3 = -4.3159552 
,4 4 = + G.4407564676X10- 3 A^" 1 
^ 5 = - 3.5663884587 X10- 6 A 6 r 2 



*-V: 



N-M-l 



(16) 



where N is the number of points and M is the 
number of constants, is 1.52 mK. The estimate of 
the standard deviation of a predicted point varies 
from 0.75 mK near the triple point, to 0.34 mK in 
the middle of the range, to 0.46 mK near the upper 
end. The values and the calculated curve are shown 
graphically in figure 5. The least squares computer 
program is given in appendix 6 and the computer 
printout of the least squares solution is given in 
appendix 7. 

10. Study of IPTS Realization Above 300 °C 

Above 200 °C, the values of temperature deter- 
mined from 3 standard PRT's differed from one 
another at the same temperature by more than might 
be expected from the calibration results. The average 
deviation from the mean of their temperature 



300 



400 






t — — i r 

H, t.p T H, b.p. 


Sn f.p. 


1 


( 


T" 

Zn f.p. 
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- 


040 
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D 


- 


060 








^Bv 


080 


1 1 . , , L 






1 



500 

68 



600 



700 



Figure 5. Deviation of International Practical Temperatures from thermodynamic temperatures 

in the range from 273 to 730 K. 

The estimated uncertainties at the defining fixed points published in the original text of the IPTS-68 are shown by the 
error bars. Key: □ high pressure range, ♦ medium pressure range, + low pressure range. 
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values versus f 68 was derived from columns 2, 3, 4, 
and 5 of tables 2, 3, and 4 and is shown in figure 6. 

AVERAGE DEVIATION OF THREE P.R.T.'s vs. t R « 




Figure 6. The average deviation of three PRT's versus the 
average IPTS temperature. 

Subsequent to all the reported gas thermometer 
measurements, we carried out special measurements 
at 408 and 414 °C to study this apparent contradic- 
tion. It has already been noted that the triple point 
resistances of the thermometers usually increased 
with time of immersion in the stirred liquid thermo- 
stat of the gas thermometer. But we have also 
observed that the triple point resistance, if already 
in substantial excess of its calibration (annealed) 
value, decreases initially when placed in the thermo- 
stat at 457 °C. We now interpret this fact as resulting 
from a balance between the rate of increase of the 
resistance from work hardening of the platinum by 
the mechanical acceleration from the stirring and 
the rate of decrease of the resistance by annealing 
of the platinum at 457 °C. At somewhat lower 
temperatures (350 to 420 °C) the rate of increase 
from work hardening remains high, but the annealing 
rate is substantially smaller, so that after the 
thermometers had been in the thermostat, a rise in 
the triple point resistance was observed. We found 
the triple point resistances to increase progressively 
with succeeding periods in the thermostat at 408 °C, 
and, although the average temperature remained 
nearly constant — both in fact and as calculated from 
the average of the values determined from the 
thermometers — the average deviation from the 
mean increased successively from 1.63 ml( to 1.89 
mK to 2.13 mK. When annealed before each measure- 
ment, the average deviation from the mean became 
0.65 mK and remained essentially constant with 
further cycles. The precision of the values of / 68 was 
the same after }{ h anneals as after 2){ h anneals. 
Although we also used different annealing temper- 
atures, the best choice, on the basis of Berry's 
isochronal step depression results [13], is probably 
450 °C. Some conclusions one may draw: 

A. Platinum resistance thermometers become 
very sensitive to mechanical shock at temperatures 
in excess of 350 °C; and up to 457 °C, the calculated 
temperature determined by them tends to drift in 



value with progressive work hardening. For rela- 
tively minimal lengths of exposure these increases 
were limited to steps of resistance corresponding 
to less than 0.30 mK. 

B. The direction of change of the resistance of a 
PRT at the zinc point from its annealed value 
because of work hardening is unpredictable on an 
a priori basis. For the thermometers we happened to 
use, the average of the values of t &8 calculated from 
the measurements were a good approximation 
(within 1 mK) of the true value of t &8 , whether 
determined from the annealed thermometers, or from 
thermometers of the same group which exhibited the 
effects of work hardening. 

C. The values of t &s calculated for a given tempera- 
ture as measured with a group of thermometers kept 
as closely as possible in an annealed condition, are in 
good agreement with one another. The best choice of 
annealing temperature is probably 450 °C, and at 
least for thermometers only slightly work hardened, 
a half-hour of annealing is sufficient. We believe 
under our conditions of measurement, the thermom- 
eters should be annealed each time before measure- 
ment is made, and, as has become our practice, the 
time during which the thermometers are subjected 
to mechanical acceleration should be kept as short 
as possible. 

Other causes of measured temperature difference, 
investigated and found insignificant, were possible 
temperature gradients in the gas thermometer bulb 
case and possible imprecision of the thermostat 
temperature. 

11. Estimation of Total Uncertainty 

The stated total uncertainty consists of the limits, 
at 99 percent confidence level, for the random errors; 
and the systematic errors, estimated conservatively 
enough to warrant about the same confidence level. 
We discriminate between errors which fall into three 
categories : 

A. Sources of random error which contribute to 
the imprecision observed in the results. 

B. Sources of error which produce a constant bias 
in the results, but which are estimable in terms 
of random errors determined from other ex- 
periments. 

C. Sources of "systematic" error which produce 
a constant bias in the results, but for which 
there is inadequate experiment and theory to 
permit their evaluation by accepted statistical 
techniques. 

The estimates of the values of the random errors 
(A and B) are expressed initially as one standard 
deviation. 

The error of the reported quantity, the deviation 
of international temperatures from thermodynamic 
temperatures, T/K-TJK &8 , can be estimated by 
combining in quadrature the errors of realizing the 
international temperatures with the errors of realiz- 
ing the thermodynamic temperatures. The errors 
for each depend upon the range of temperature, and 
for the thermodynamic temperatures, also upon the 
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range of pressures. There are two values of tempera- 
ture for which a larger number of values was de- 
termined — 373 and 730 K. We propose to assess 
the errors at these two temperatures, being values 
at or near the extremes of the range, as well as to 
consider the statistical evaluation of the imprecision 
derived from fitting the curve in terms of the esti- 
mate of the standard deviation of a predicted point. 
The errors in the realization of the international 
temperatures are those from: 

(1) The calibrations. The systems of measuring 
instruments, fixed points and thermometers is 
sufficiently exact that for the four standard 
PRT's used in the measurements, the con- 
stants A and B in the equation 



R=R (1+At'+Bt' 2 ), 



(17) 



(2) 



determined in four successive calibrations 
over a period of two years defined the value 
of the temperature at the zinc point with an 
estimated standard deviation of the mean of 
the thermometers of 0.20 mK. 
The values of t m determined by three standard 
instruments which had the precision of cali- 
bration stated in (1) above. 
The estimate of the standard deviation of the mean 
in realizing the international temperature can be 
evaluated by combination of the estimates of the 
standard deviations of the mean of the temperature 
imprecision of the 3 thermometers because of the 
imprecision of the constants of calibration, Si; and 
3 factors involved in the measurement of a tempera- 
ture: the imprecision of the thermostat bath regula- 
tion, S 2 ; the imprecision of the measurement of the 
resistances, S 3 , and the uncertainties because of the 
effects of work hardening discussed in section 10, £ 4 . 
These are given in table 7. 

Table 7. Errors of realization of Us estimated standard 
deviation of the mean 









Type of error 


Element 


^68 
°C 




















A 


B 


O 






DiKfis 


mKfls 


mK 68 


Si (Calib.)-- 


457 
100 




0. 45 
.30 




S 2 (Thermo.) _ 


457 
100 


0. 20 
. 10 






aS 3 (Resis.)--- 


457 


. 02 




0. 10 b 




100 


. 02 




. 05 


Si (Hard.)-— 


457 
100 






. 33 

. 07 


V^s? 


457 


. 21 


. 45 


. 34 




100 


. 12 


.30 


.086 



a These values are not to be added with A and B on the 
basis that their combination with random errors is not philo- 
sophically acceptable. They are stated at 1/3 the value we 
estimate, for easy comparison with the errors of Type A and 

b Thesc are values, suggested by experiment, for the ac 
effects of shunting by dc leads, i.e., a lossy capacitor. 



The approximate gas thermometer temperature 
was calculated from eqs (2) through (5a) and the 
thermodynamic temperature was then found by ac- 
counting for the effects of thermomolecular pressure 
and gas imperfection. The approximate gas ther- 
mometer temperature can be expressed in terms of 
the ratio in a way useful for error analysis as 

T f GT /K / G T=T l /K(Z 2 /Z l ) 



(Ti/K)T 2 



£*(1-f!©' v» 



The value of the second term is small relative to 
T GT , so that the relative errors of the products in 
the first term are directly reflected as relative errors 
in the gas thermometer temperature. In the order of 
the quantities in eq (18), the errors are 

(1) The relative uncertainty of the pressure ratio, 



mm- 



The error has two main components, an imprecision 
depending upon the regulation of the pressure and 
the variability of the diaphragm, and a fixed uncer- 
tainty that can be evaluated from the imprecisions 
in p, g and It, The pressure difference between the gas 
thermometer and the manometer was measured by 
the diaphragm transducer before and after the bath 
temperature readings. The estimates of the standard 
deviation of the mean were calculated from the devia- 
tions for the experimental pressures at 100 and 
457 °C and also at the fiducial point. 

For 457 °C, the estimate of the standard deviation 
of the gas thermometer temperature from this source, 
for a pressure of 10 5 Pa is ±0.23 mK GT . To this 
should be added a possible random fluctuation of the 
regulated pressure, estimated in terms of a variance 
from the manometer setting. This is estimated to 
amount to ±0.033 Pa at 1 at in (101,325 Pa) and 
0.02 Pa at 3.5X10 4 Pa, These effects combined in 
quadrature with the variability of the diaphragm 
are, for one standard deviation, 0.33 mK GT at 10 5 Pa 
and 0.52 mK GT at 3.5 X10 4 Pa. 

For 100 °C, the estimate of the standard deviation 
of the mean of the diaphragm readings affects the gas 
thermometer temperature by only 0.34 mK GT for 
10 5 Pa, and 0.17 mK GT at 1.8 X10 4 Pa. An estimate 
of random fluctuation in gas thermometer tempera- 
ture, from the manometer setting at 10 5 Pa of 
0.033 Pa is equivalent to 0.12 mK GT and at 1.8X10 4 
Pa of 0.02 Pa is equivalent to 0.41 mK GT . These 
errors combined in quadrature are 0.12 mK GT for 
10 5 Pa and 0.44 mK GT for 1.8X10 4 Pa, 

The total uncertainty of the pressure ratio at the 
manometer is discussed in [5]. It was stated to be 
±1.5 ppm at the 99 percent confidence level, or 
about 0.5 ppm for 1 standard deviation of the 
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220-208 0-77-2 



mean. The error can be evaluated from measured 
quantities, so is of a "B" type. 

(2) The contribution of the pressure head to the 
measured pressure in the present temperature range 
varies between 4.5 ppm and 11.5 ppm and is un- 
likely to be in error by more than 2 or 3 percent. 
Therefore, the error of the pressure head correction 
is relatively insignificant. 

(3) The errors in the thermal expansion measure- 
ments used to calculate the thermal expansion of 
the bulb. It is believed that the systematic effects 
are relatively small. Therefore, the error can be 
evaluated from the imprecision of the measurements, 
and is of type "B". The volume expansion is calcu- 
lated from 



where 



v/Vo=(i+w Aty 

-i 1 l~k 
a = T -—-> 



(19) 
(20) 



l is the length at °C, and t Q8 is the temperature 
at which I was determined. As stated in section 6, 
a fifth order power series was fitted by least squares 
to the data. The estimate of the standard deviation 
of a predicted point varied from a minimum of 0.8 
part in 10 7 from to 54 °C, 1.9 part in 10 7 at 457 °C, 
and was 1 part in 10 7 at 100 °C. Therefore, the 
error contributed to the thermodynamic tempera- 
ture, expressed as the estimate of a standard devia- 
tion is 0.38 ppm at 100 °C and 0.62 ppm at 457 °C. 
(4) The uncertainty in the dead space calculation 
is primarily due to the uncertainty in the determina- 
tion of temperatures along the connecting tube. 
These were measured b} r thermocouples, the measure- 
ment problems of which have been discussed earlier 
in this paper. The effect of the temperature uncer- 
tainty can be evaluated from an approximation for 
the value of Vt, which is 



Vi 



Yi 



=S*^ccV/T 



(21) 



where T=(*+23)/2+273.15 is the absolute tem- 
perature averaged between the temperatures of the 
room and the thermostat. We have evaluated the 
error in the gas thermometer temperatures resulting 
from a total uncertainty in Tof 8t= ±0A%X(t— 23) 
^C. For instance at 457 °C, the average temperature, 
T, is assumed uncertain by 1.74 °C for the measure- 
ments (corresponding to an uncertainty of about 
2.5 °C at t itself). A summary of the error calcula- 
tions, which become insignificant below 141 °C, is 
given in table 8. 

The value of one standard deviation of the mean, 
S m , is about % of the total uncertainty given in the 
final column. 

The remaining sources of error come from un- 
certainties in the evaluation of thermomolecular 
pressure and gas imperfection effects. The thermo- 
molecular pressure correction was largest for states 
10.01 and 10.02 at 457 °C, where it amounted to 



16.4 mK GT . The only significant error in the correc- 
tion derives from the imprecision of the measure- 
ments for determining thermomolecular pressure. 
The estimate of the standard deviation of the mean 
for the results at each temperature measurement 
was 1.1 percent at 140.5 °C, 0.68 percent at 260 °C, 
1.1 percent at 374 °C and 0.92 percent at 457 °C. 
The error is of type B, and amounts to about 1 
percent of the calculated correction, which in general 
is not significant. 

The effects of gas imperfection were calculated 
from the virial coefficients for helium, selected as 
we described earlier. The largest correction for gas 
imperfection was made for states 9.04 and 9.05 
at 400 °C, but it was almost as large for the numerous 
measurements at 457 °C with ^ = 10 5 Pa. We 
postulate that the instrument error in the measure- 
ments to evaluate Z=pV/ET=l J rBp/RT produces 
a total uncertainty in Z of at least ±5 ppm, and is 
that small only for exceptional work. On the basis 
of a constant error in the virial coefficient of ±0.11 
cm 3 /mol, this in turn makes an error for the com- 
bination of B 2 — Bi of ±0.16 cm 3 /mol, which is 
perhaps that small only in the range 0-150 °C. At 
least twice that total uncertainty should be assumed 
for the measurements at 400 °C (states 9.04 and 
9.05), and for the measurements at 457 °C. These 
uncertainties amount to ±4.4 mK for states 9.04 
and 9.05 and ±3.7 mK for the higher pressure 
measurements at 457 °C. 

We can summarize the errors in the thermo- 
dynamic temperature as being 

S m (T)=S m (GT)+S m (TheTmp)+S m (Imp) (22) 

where S m (GT) is the estimate of the standard 
deviation of the mean of the gas thermometer 
temperature, S m (Thermp) is the estimate of the 
standard deviation of the mean of the thermo- 
molecular pressure correction, and S m (Imp) is 
the estimate of the standard deviation of the mean 
of the gas imperfection correction. 

The estimate of the standard deviation of the 
mean of the gas thermometer temperature can in 
turn be broken up into 

S w (GT)=#„/PR)+S m (TE)+S m (DS) (23) 

where S m (PR) is the estimate of the standard 
deviation of the mean of the gas thermometer 
temperature produced by errors in the pressure 
ratio. 7 

S m (TE) is the estimate of the standard deviation 
of the mean of the gas thermometer temperature 
produced b}^ errors in the values of the thermal 
expansion of the bulb. 

S m (DS) is the estimate of the standard deviation 
of the mean of the gas thermometer temperature 
produced by errors in the evaluation of the cor- 
rection for the deadspace. 

The total can be found b}' addition in quadrature 
of the values given in table 9. 

7 The error of the pressure head is estimated to be insignificant. 
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Table 8. — Estimate of the uncertainty, 8toT, in T GT from the assumed uncertainty, 
8t, in Thermocouple readings 



t 


8t 


T 


8T/T 


T l T 2j/ P 2 


8 Tgt 


°Ce8 


K(, s 


Kq$ 




Kqt 


Kgt 


457 


1. 74 


515 


0. 0034 


0. 7583 


0. 0026 


425 


1.61 


499 


. 0032 


. 7124 


. 0023 


392 


1.48 


482 


. 0031 


. 6670 


. 0021 


360 


1.35 


466 


. 0029 


. 6248 


. 0018 


327 


1. 22 


450 


. 0027 


.5864 


. 0016 


295 


1. 09 


434 


. 0025 


. 5438 


. 0014 


263 


.96 


418 


. 0023 


. 5027 


. 0012 


232 


.83 


402 


. 0021 


.4621 


. 0010 


200 


. 71 


386 


. 0018 


. 4224 


. 0008 


173 


. 60 


373 


. 0016 


.3995 


. 0006 


141 


.48 


357 


. 0013 


. 3522 


. 0005 



Table 9. Sta7idard deviations of the mean, and systematic errors in the realization of thermodynamic temperatures 







Class of error 




A (mK) 


B (mK) 


O (mK) 


>S' m (PR)_ 


* ' °\3.5X10 4 Pa 

/1 05 Pa 
100 °C\1.8X10 4 Pa 

457 °C 
100 °c 

457 °C 

100 °c 

457 o C /10»Pa 

40 ' U \3.5X10 4 Pa 

100 o r /10* Pa 

lUU U \1.8X10 4 Pa 

f 4 5 7 o C flO*Pa 
140/ ^| 3 .r )X 10 4 Pa 

ioo°c/ 105Pa 


0.33 
. 52 
. 12 
. 44 

Negligible 

. 33 
. 52 
. 12 
. 44 


0. 37 

. 18 

. 45 
. 14 

. 58 
. 23 




£ m (TE) 




S m (DB) 

£ m (Thermp) 

£ TO (Imp) 

Vz,-*s* 


0. 87 
. 10 

1. 23 
. 43 

. 65 
. 12 

1. 51 

. 97 
. 66 
. 16 



ft The errors of Classes A and B are given as one standard deviation. The estimated values of Class 
C are divided by 3 in order to give a "common basis" of comparison with A and B. 

Thus the estimate of the standard deviation of The observed imprecision should be consistent 
the mean of T/K—Tq 8 /K 6S is the combination of with the sum of "A" of table 10 and S 4 of table 7 
A, B, and C from tables 7 and 9, as follows: (the observed systematic deviation from work 

Table 10. Estimate of the standard deviation of the mean and of systematic errors in the determination of T/K-Tm/Km 







Class of error 




A (mK) 


B (mK) 


C (mK) 


457 °C 

100 °c 


10 5 Pa 

3.5 x 10 4 Pa 

10 5 Pa 

1.8 x 10 4 Pa 


0.391 
.56/ 

.171 
.46/ 


0. 73 
.38 


1. 55 
1. 03 

. 67 

. 18 
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hardening of the PRT's). These latter are, in quadra- 
ture, the following: 

5.2 X10" 4 at 10 5 Pa and 457 °C 

6.6X10" 4 at 3.5X10 4 Pa and 457 °C 

1.9 X10" 4 at 10 5 Pa and 100 °C 

4.8X10" 4 at 1.8X10 4 Pa and 100 °C. 
The corresponding "observed imprecision", the 
standard deviation of the mean of the experimental 
values of T/K—T Q8 /K &8 , at each temperature is 

6.6X10" 4 for 10 5 Pa and 457 °C 

1.02X10" 3 for 10 4 Pa and 100 °C. 
The large experimental value at 100 °C may indi- 
cate sources of random error that were not adequate- 
ly appraised (as discussed in [8]). The observed 
imprecisions can also be evaluated more broadly 
in terms of the deviations of the 72 values of 
T/K— T 68 /K 6S from the least-squares fit. The stand- 
ard deviation of a predicted point for T/K—T 6S IKq$ 
is4.6Xir 4 at457°C and 3.5X10" 4 at 100 °C. 

The estimates of the total uncertainty can be 
derived by adding the estimates of type A and B in 
quadrature, and multiplying by 3 as the approxi- 
mate factor from Student's Table for a 99 percent 
confidence level. The systematic errors, part C, 
are similarly increased, because they were arbitrarily 
stated at % of their estimated values in order to 
compare them with the values of Types A and B. 
The summed values for 457 and 100 °C are given 
in table 11. 



Table 11. 



Total uncertainty of T/K—T^/K^ at 730 and 
373.15 Km • 



457 °C 


10 5 Pa 


2.7X 10~ 3 random error 
4.6X10 -3 systematic error 




3.5X10 4 Pa 


2.8 X 10~ 3 random error 
3.1X10 -3 systematic error 




10 5 Pa 


1.2X10 -3 random error 


100 °C 




2.0X10 -3 systematic error 




1.8X10 4 Pa 


1.8X10 -3 random error 
.54X10 -3 systematic error 



12. Discussion 

All data, including those published previously, 
have been interpreted in the light of further informa- 
tion and improved understanding of the apparatus. 
The revaluation of the volume of the bulb makes an 
appreciable difference in the dead space calculation 
for the data published previously (amounting to as 
much as 2.8 ml( at 414.75 K). In contrast, the 
difference in composition of the top of the bulb from 
the remainder causes no net change in bulb volume 
with temperature (compared with an all 88 percent 
Pt-12 percent Rh composition) in excess of 1 ppm, 
because the shape of the distortion is compensating 
in the volume to the first order. The problem is 
minimized because the total thermal expansion of the 
two alloys becomes equal for £~360 °C. Inasmuch 
as we had used the values for 88 percent Pt-12 
percent Rh in the previous calculations there is no 



error exceeding 1 ppm from this cause in any of the 
results, old or new. 

All the evidence of system cleanliness — RGA 
scans, diaphragm, and gas thermometer stability, 
reproducibility before and after gas thermometer 
bakeouts — supports the proposition that the results 
are free of significant effects of sorption (except set 
G which was not included in the final results). 
Above 0.01 °C, gas thermometer temperatures are 
increased by the effects of sorption. It is not generally 
recognized (although it was mentioned in [9]), that 
our earliest measurements yielded " thermodynamic 
temperatures" in excess of 373.15 K at the normal 
boiling point of water; these values declined as the 
gas thermometer became cleaner. We believe we 
could adjust the purity of the gas thermometer to 
produce a selected value of the steam point over a 
range of about 0.1 K. In the less clean condition, we 
expect that there would be less precision and stability, 
but the gas thermometer would probably be decep- 
tively stable even with enough sorption effect to 
give"T th "=373.15K. 

Operating in the higher temperature range em- 
phasized certain problems, viz., the measurement 
and treatment of dead space temperatures and of 
thermomolecular pressures, the evaluation of the 
effects of gas imperfection, and the realization of 
international temperatures, all of which can be 
handled better in the future with higher resulting 
accuracy. We have installed a new gas thermometer 
that has a much improved thermocouple installation 
for measuring the dead space temperatures. Special 
annealing and testing equipment was built to assure 
that these couples were uniform and had emfs 
in close agreement with accepted values. The effects 
of thermomolecular pressure and gas imperfection 
are both encountered in any of our gas thermometer 
measurements, but they can be evaluated with little 
loss of accuracy by a modification of the customary 
measuring techniques : 

(1) Independent measurement of thermomolec- 
ular pressures. The results of some measurements 
have been reported. Further work is planned where 
conditions are more exactly defined and more ac- 
curate values are obtained at lower average pressures. 

(2) Use of relatively small fiducial pressures. Our 
ability to determine the virial coefficients is limited 
by our ability to determine thermomolecular pres- 
sures. Therefore, the highest accuracy results when 
only one factor is involved. The measurement of 
thermomolecular pressures and the theoretical justi- 
fication for the transfer of results to the gas thermom- 
eter will be the subject of a separate paper. 

It has recently become possible to operate the 
gas thermometer with no loss of accuracy in a pres- 
sure range much lower than was the limit ten years 
ago. Of special benefit in this respect is the ten-fold 
improvement in the calibration of gage blocks. No 
new manometer errors become significant in a low 
range of pressure based on a fiducial value of 2.54 cm 
of Hg. The great advantage of this lower pressure 
range is that the correction for gas imperfection is 
sufficiently small that the uncertainty is less than 
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1 ppm at 457 °C for a pressure of 6.78 cm Hg, and 
the thermomolecular pressure correction, though 
large, is measurable with such accuracy that its 
uncertainty is about the same. 



13. Summary and Conclusions 

The deviation of international temperatures from 
thermodynamic temperatures has been evaluated 
by a constant volume helium gas thermometer be- 
tween and 457 °C. Previously reported values were 
modified for a correction in the dead space effect, 
and also in light of a reinterpretation of the effects 
of gas imperfection and thermomolecular pressure. 
The results are consistent with independent measure- 
ments of each of these effects. As examples, the total 
uncertainty, consisting of the estimate of the ran- 
dom errors at the 99 percent confidence limits and 
the systematic errors, intended to be conservatively 
enough stated to warrant about the same confidence, 
was calculated for each of two pressures at 100 and 
457 °C. The estimate of the systematic error in each 
case was appreciably larger for the higher pressure 
range, but the results agree with those measured in 
a lower pressure range. Hence the uncertainties for 
the lower pressures are regarded as the proper ones, 
and are given following the value of the deviation: 



t 



457 °C 
100 °C 



T/K — T^/K^ 



-0. 0794 ±2. 8 X 10- 3 random error 

±3. 1 X 10~ 3 systematic error, and 

-. 0252 ± 1. 8 X 10- 3 random error 

±. 54 X 10~ 3 systematic error 



The errors at temperatures between 100 and 457 °C 
ought to be approximately proportional to the 
temperature so that one obtains : 

At the tin point, * fl8 = 231.9681 °C 



T/K-T 9f JK w 



-0. 0439 



±2.2X10 -3 random error 
± 1.5X 10 -3 systematic error 



At the zinc point, * fl8 = 419.58 °C 



T/K — T 6 sIKqs 



-0. 0058 



±2.8X10~ 3 random error 
±2.8X10 -3 systematic error 



Improved equipment and improved procedures 
are expected to substantially lower the uncertainties 
in the next set of measurements, which have been 
commenced. 
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APPENDIX I 

100 * FORTRAN IV PROGRAM FOP CALCULATION OF IPTS-68 VALUES 

110 * FROM PRT RESISTANCE MEASUREMENTS* NAMED "BRIDGE" 

120 DOUBLE PRECISION FTC 5>« PC 5> « V« AC 10)«BC 10) * TC 5)« TIjZ 

130 DIMENSION NC 5) 

140 CALL DEFINE C 1., 7HTCAL74.. ) 

150 CALL DEFINE <2*4HPFT, ) 

160 READ C 1* 1 10)CACJ).,BCJ),J = 1, 10) 

170 110 F0RMATCD15.8>2X,D15.8> 

180 120 READ C2« 130>CS1*S2jCNCJ>«FTCJ>*FCJ>«Js1«4>> 

190 130 F0FMATC2X,F6.2.,X.,F6.3.»C4C/2X>I4.,D15.8>D15.8))) 

200 WRITE C9> 130)CS1> S2,CNCJ).,FTCJ).,FCJ),J=l/4) > 

210 WRITE (9* 150) S1/S2 

220 150 FORMATC/UX* ' STATE'F6. 2.» 2X* F6« 3* /) 

230 DO 240 J=l,4 

240 W=PTCJ)/CPCJ)-0.996D-03) 

250 I=NCJ)-700 

260 TCJ)=C-ACI)+DSQFTCACI )**2-4.* C 1 . -V)*BC I ) ) ) /C 2.* EC I ) ) 

270 Z=0.45D-03*TCJ)*CTCJ)/100.- 1 . )*C TC J ) /4 19 . 58- 1 . )*C TC J) / 630. 74- 1.) 

280 TCJ)=TCJ)+Z 

290 WRITE C9*230)MCJ)*TCJ) 

300 230 FOFMATC 'THERM. NO. ' I 3* 3X> ' TC I NT. 1968) =' D17.10) 

310 240 CONTINUE 

320 Tl a CTC2) + TC3) + TC4))/3 

330 WRITE C9* 270) Tl 

340 270 FORMATC 'AVE. TCINT. 1968) = ' D 17 . 1 0* /» /) 

350 GO TO 120 

360 END 

TCAL7U 

0.39849575D-02 -0. 587624 1 5D-06 

0.39850801D-02 -0. 58772 1 19D-06 

0.39848143D-02 -0. 53762 18 5D-06 

0. 398499 17D-02 -0. 58769 654D-06 

0. 39849663D-02 -0. 58756733D-06 

0.39847884D-02 -0. 58757724D-06 

0.39848898D-02 -0. 58760947D-06 

0.39850236D-02 -0. 58778 636D-06 

0.39838174D-02 -0. 58757363D-06 

0.39855387D-02 -0. 58755669 D-06 

PRT 

10 9.01 10.036 

12 710 0.27587830D 02 0.25557448D 02 

14 701 0.63945402D 02 0.25577651D 02 

16 703 0.63731295D 02 0.25492645D 02 

18 707 0.63956557D 02 0.25582540D 02 
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100 PEM PYE* IN XBASIC* 9.29*75 

110 PEM USING TWO FILES* DATA AND LSFIT1 

1 1 5 DBL 

120 READ Cli C2*C5*C6 

130 DATA 8.74E-6* 3E-9* 13.0105E-6* .384333E-08 

140 DIM V C 18)* V<6j 18)* G( 13)* E( IS) 

15 0PEADB0*B1*B2*B3*B4*NC 1)*NC2)*NC3) 

160 DATA 135* . 54 58 46E-02* . 1 1 349 7E- 04* - . 152447E-0 7* .9 06 03 3 E- 1 1 

17 0DATA1*4* 12 

1S0FOPI=1TO13 

190PEADGCI ) 

200NEXTI 

210 DATA 62.2*89.6* 1 16-9* 144.2* 17 1*6* 198*8* 226.3 

220 DATA 253.6*280.9*308.3*335.6*363.0*390.3 
230 FILES DATA**DISP 

240 PEAD #1* S0* S1*T9*E0* El* E2 
242 IF END #1 THEN 109 

250 PPINT* 260* S0* S 1* T9*E0* E1*E2 

260 FMT F7.3*F7. 2*F1 1.5* 16* I 5* I 5 

270 FOP 1= 1TO 13 

28 PEAD #1*E(I) 

290 PPINT* 300* EC I )* 

300 FMT 15 

310NEXTI 

320 PPIN^ 

330 D2=0 

34 V=0 

35 V0=B1*T9 + B2*T9 t 2+B3*T9 t 3+B4*T9 t 4 
360 FOP I=1T013 

370 T1=B0 

38 ECI)=1 13+CE0-787)*. 147+ EC I >*1. 0089 5 

390 ECI)=ECI )* IE- 3 

400 LET F=B1*T1-EC I )+B2*Tlt 2+B3*T 1 1 3+B4*T 1 14 

410 LET G=Bl + 2*B2*Tl + 3*B3*Tlt2+4*B4*Tlt3 

420 T1=T1-F/G 

430 IF ABSCF/G)>2E-4THEN400 

440LETD1=T1+CT1/100)*C-. 006+ CT1/ 100- 1)*C . 04 1 06- 7 . 363E- 5*T 1 ) ) 

450LETYCI ) = 1/CD 1+273. 15) 

460NEXTI 

470 CALL LSFIT1 

480 FOP 1=62 TO 144 

490 LET X=BC 1)+BC2)*I+BC 3)*I t2 

500 GOSUB 930 

510 NEXT I 

520 FOP 1= 145T0363 

530 LET X=BC4)/CIt 2) +BC 5) /I + BC 6) +BC 7) *I + BC 8 ) *I t2 

54 GO SUB 9 30 

550 NEXT I 

560 FOR 1 = 364 TO 390 
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570 LET X=Y< 12)+CY< 13)-Y( 1 2) )*( I -363) /27 

580 GOSUB 930 

59 NEXT I 

600 IF V0-E<13)>.1 THEN 700 

610 FOR 1 = 39 1 TO 472 

620 LET X=Y< 13)+( l/(T9 + 273. 15)-Y< 1 3) )*( I -390) /82 

630 GOSUB 9 30 

640 NEXT I 

650 FOP 1=473 TO 688 

660 LET X= l/(T9+273. 15) 

670 GOSUB 9 30 

68 NEXT I 

69 GO TO 78 

700 FOP 1 = 39 1 TO 526 

710 LET X=Y< 13) + < l/(T9 + 273. 15)-YC 1 3) )* ( I -390) / 136 

720 GOSUB 9 30 

730 NEXT I 

740 FOP 1=527 TO 688 

750 LET X= l/(T9 + 273. 15) 

760 GOSUB 930 

770 NEXT I 

780 LET D2=D2-55*(C5-C1+(C6-C2)*T3)*T2 

790 T( l)=<E0+El-787)/40+20 

800 T(2) = (E0+E2-787)/40+20 

810 FOR I=1T02 

820 T1=T(I) 

8 30 LET Dl=Tl+(Tl/100)*(-0.006+(Tl/100-l)*(0.04106-7.363E-5*Tl)) 
840 T<I)=Dl+273. 15 

850 NEXT I 

860 V-W+ 13/T(2) + 2*(42-D2)/(T<2) + T( l)) + 30/<T( 1)+1/Y( 1)) 
870 LET F=4*980. 6/(82. 07*1013250) 
880 PRINT 

890 PRINT" BOOK STATE PI DI SPL" 

900 PRINT S0; Sl;F*V*.09 29*D2*.9 29 
910 PRINT 
911 WRITE #3*D2*.9 29 
920 GO TO 240 
930 LET T2=l/X-293. 15 
940 LET T3=T2+20 

950 LET W=V+X*( 1 + (C1 + C2*T3)*T2) 
960 LET D2=D2+(C5-C1+(C6-C2)*T3)*T2 

9 70 RETURN 

980 DATA 10.036*9. 05* 400. 07842.. 1041* 50* 14 

990 DATA 79*207* 488*849* 1266* 1721*2137*2458* 2671 

1000 DATA 2846*2979*3037*3073 

109 STOP 

1100 END 



726 



APPENDIX III 



460LETA9=1 
470FOPA=2TO3 

48 0FOPI=N(A- l)TON(A) 

49 01 FN < A) -N ( A- 1) = 3TH EN 520 
500LETM=5 

510GOTO530 

520LETM=3 

530FOPK= 1TOM 

540iFM = 3THEN 57 

550LETVCK>I)»G(I>t(K-3) 

560GOTO 580 

570LETV(K*I)=G(I)t(K-l) 

580NEXTK 

59 0NEXTI 

600LETZ = N<A> 

610GOFUB68 

620FOPL=1TOM 

630LETB( A9)=-V(M+ 1.,7,+L) 

640LETA9=A9 + 1 

650NEXTL 

660NEX"TA 
670 RETURN 

68 0FOPI=N(A- l)TON(A) 

690LETV<M+ 1*I)=Y(I ) 

700NEXTI 

710FOPI=N<A>+1TON<A>+M 

720FOPK=1"-OM+1 

730LETV<K,I)=0 

740NEXTK 

750LETV(I-Z>I)=1 

760NEXTI 
770FOPK= 1TOM + 1 
78 01 FK=M+ 1 THEN 8 10 

79 0LETT=1 

800GOTO840 

8 10LETT=- I 

820GOTO850 

830LETT=-T 

840IFK=1THEN910 

8 50FOPL= 1TOK-1 

860LET?(L)=0 

8 70FOPJ=N(A- l)TCN(A) 
880LETS<L)*S(L) + V<K*J)*V<L.»J) 
89 0NEXTJ 

9 00NEXTL 

9 10FOPI=NCA-1)TON(A)+M 
9 20IFK=1THEN 9 60 
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930FOFL=1TOK-1 

940LETV<K*I)«l'<K*I>-S<L>*V(L*I > 

9 50NEXTL 

9 60MEXTI 

9 70IFT<0THEN 1070 

980LETD<K)=0 

99 0FOPI «N< A- 1 ) T0N( A) 
1000LETD<K)=D(K>+V<K* I >*IKK* I > 
1010MEXTI 

1020LE" P D=SCP(D(K>> 
1030FOPI=N(A-1)TON(A)+M 
1040LET\/(K*I>=V<K*I)/D 
1050NEXTI 

10601 FT>0TH EM 8 30 
1070NEXTK 
108 0PETURN 



DATA 



10 10.018* 10.08*4 57.27 365*647*54* 

12 86*223*51 1*9 15* 1377* 19 00* 239 2* 277 2* 3017*3243*3421*3548*3665 

14 10.018* 10.09*4 57.27530*647*54* 

16 8 3*219*507*9 16* 1376* 19 03*2401* 278 7*3033*32 57*3431*3545*3654 

18 10.018* 10. 10*457.274 17*639*54* 

20 9 0*230*519*9 24* 1376* 189 0*237 3*2753*2998*3228*3410*3531*3649 

22 10.018* 10. 1 1*457.27320*639*54* 

24 88*227*512*9 19* 137 1* 189 4*2387*2772*3019*3244*3421*3535*3646 

26 10.018*9.06*424.9 6377*627*49*3 

28 78*200*455*8 29* 1227* 1698*2145* 2498*2733*2944*3 115*3221*3314 

30 10.018*8.07*39 1.74 159*587*45*7 

32 71* 180*4 05*728* 1078* 149 1* 1889* 2209*2433*263 5*28 05*289 7*299 

34 10.018*7.04*359.57478*579*40* 10 

36 63* 160*359*640*948* 1312* 1662* 19 50* 2161*2350*25 13*2607*268 4 
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100 PEM DEDSPS, IN XBASI C, 9 . 29 .75 

110 REM WITH THREE FILES: "DATA, " GIVING S0, S 1, T9, E0, E 1, E2 AND 

120 REM Al THROUGH El; "DATA2, " LISTING VALUES OF "A" FOP I.D. 

130 REM OF TUBE, AND LSFIT1 FOP LEAST SQUARES. 

140 DBL 

150 READ C1,C2, C3, C4 

160 DATA 8.7252E-6,2.5279E-9,3.00786E-12,2.67373E-15 

170 DIM Y< 18),VC6, 18), G< 18), HC 18), E( 18),P(390) 

18 0READB0, E 1, B2, B3, B/i, N( 1 ) , N< 2) , N( 3) 

190 DATA 135, . 545846E-02, . 1 13497E-04,-. 152447E-07, .906033E-1 1 

200 DATA 1,4, 12 

210 FOR I=1T013 

220 READ H<I) 

230NEXTI 

240 DATA 62.2,89.6,116.9,144.2,171.6,198.8,226.3 

250 DATA 253.6,280.9,308.3,335.6,363,390.3 

26 FILES DATA, DATA2, DISP 

270 FOP I=1TO390 

28 READ #2, A 

290 PCI)=0. 90123+2. 748E-4*A 

300NEXTI 

310 READ #1,S0,S1,T9,E0,E1,E2 

320 IF END #1 THEN 1210 

330 READ #3, £2 

340 PRINT Z2 

350 PPINT, 360, S0, S1,T9,E0, E1,E2 

360 FMT F7.3,F7.2,F1 1.5,16,15,15 

370FORI=1TO13 

380 GU)=H(I)-Z2 

39 0NEXTI 

400 H2=0 

410 IF Z2<.5 THEN 450 

420 H2=H2+1 

430 Z2=Z2-1 

440 GOTO410 

450 PRINT 

460 FOP I=1T013 

470 PEAD #1,E<I) 

480 PRINT,490,E<I), 

49 FMT 15 

500NEXTI 

510 PRINT 

520 V0=B1*T9 + B2*T9 t2+B3*T9 r3>-B4*T9t 4 

530FORI= 1T013 

540 E(I)=( 113+CE0-787)*. 147+ E( I)* 1 . 00895) / 1000 

550 H1 = E(I) 
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560 GO SUB 1140 

570 Y(I)=l/(Dl+273. 15) 

58 0NEXTI 

590 CALL LSFIT1 

600 Dl=(E0-787)/40+20 

610 D3=E2/40 

620 D4=El/40 

630 T6=D1+D3 

640 D2=T6-23 

650 X7=( l+CCl + C2*T6+C3*T6t2+C4*T6t3)*D2)t3 

660 V5=5.367*X7/(T6+273. 15) 

670 FOP I=1T07 

68 GO SUB 1100 

69 0NEXTI 

700 FOP I=8T042-H2 

710 T6=Dl+D4+<D3-D4)*<42-H2-I)/(42-H2-8) 

720 GOSUB 1 100 

730NEXTI 

740 FOP I=43-H2T062-H2 

75 T6=l/Y( D + CD1 + D4+273. 15- 1/Y( 1 ) )* ( 62-H2-I ) /20-273. 15 

760 GOSUB 1 100 

770NEXTI 

780PRINT 

790 FOP I=63-H2T0 144-H2 

8 00 T6=1/(B( 1) + B(2)*I+B(3)*(It2))-27 3. 15 

8 10 GOSUB 1 100 

S20NEXTI 

830 FOP I=145-H2T0363-H2 

840 T6=l/(B(4)/(It2)+E(5)/I + B(6) + B(7)*I + B(8)*(It2))-273. 15 

8 50 GOSUB 1100 
860NEXTI 

870 FOP I = 364-H2TO390-H2 

88 T6=l/< (Y< 13) -Y( 12) )*< I - ( 363-H2) ) /27+YC 12)) -2 73. 15 

89 GOSUB 1100 
900NEXTI 

9 10 I = 39 

920 IF l'0-E(13)>.l THEN 1020 
930 FOP J=39 1-H2 TO 472 

940 T6=1/(Y< 13) + ( 1/CT9 + 273. 15)-Y< 1 3) )* ( J-C 39 0-H2) ) /( 82+H2) ) -273. 15 
9 50 GOSUB 1 100 
960NEXTJ 

970 FOP J=473 TO 526 
980 T6=T9 
99 GOSUB 1100 
1000NEXTJ 
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1010 GO TO 1060 

1020 FOR J = 39 1-H2 TO 526 

1030 T6= 1/CYC 13>+C l/(T9 + 273. 15)-YC 1 3) )* ( J-< 39 0-H2) ) /( 1 3 6+H2) ) -273 . 15 

1040 GO SUB 1100 

1050NEXTJ 

1060 PPINT" BOOK STATE SUM V*TAU" 

1070 PPINT S0;Si;V5 

1080 PPINT 

109 GO TO 310 

1100 D2=T6-23 

1110 X7=( l+(Cl + C2*T6+C3*T6t2+C4*T6t 3>*D2> t3 

1120 V5=V5+.78 539 5*P(I ) t 2*. 9 29*X7/ < T6+273 . 15) 

1130PETT7PN 

1140 T1=B0 

1150 F=E1*T1-H1 + B2*T1 t 2+B3*T 1 f 3+B4*T 1 1 4 

1160 LET G=B1 + 2*B2*T1 + 3*B3*T1 t2+4*E4*Tl *3 

1170 T1=T1-F/G 

1180 IF AES(F/G)>2E-4 THEM 1150 

119 Dl = TH-CTl/100)*(-.006+CTl/ 100- 1 ) *< 0. 04 1 06- 7 . 363E- 5*T 1 ) ) 

1200PETUPN 

1210 STOP 

1220 END 
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100 * PROGPM "S-UMMA"j REVISED APRIL 21*1975 

110 * FORTRAN IV DOUBLE PRECISION, BULB I 1 1 B, DIAPHRAGM I A 
120 DOUBLE PRECISION C* C 1* C2, C3# C4# V* T* El, E2, E3, E4 
130 DOUBLE PRECISION H,H 1, T9, PI , V9, D U C5, C6 
140 DOUBLE PRECISION T 1, T2, T3, T4> P, Y, Z, V 
150 DATA C, V/0. 17255D-03, 0.4325D 06/ 
160 DATA Cl,C2/0.8750894565D-05,0.2278241267D-08/ 
170 DATA C3,C4/-0.2476697519D-1 1, 0. 346 19 03664D- 14/ 
175 DATA C5/-0. 19 1 1694941D-17/ 
180 CALL DEFINE <1,6HDATA3,) 

190 60 READC 1,65) S0, S 1, El, E2, E3, E4 
200 65 F0RMATC4X,F6.3,X, F6.2,X, D9.2,X, D9 . 2,X, D9 . 2,X, D9 . 2) 

210 WRITE C9,65)S0, S 1, E 1, E2, E3, E4 

220 READ ( 1,75) T,H,H 1 

230 75 FORMAT ( 4X, Dl 3* 6,X, D 1 3. 6,X, D 10. 3) 

240 WRITE <9,75)T,H,H1 

250 READ (1,85) T9,PI,V9,D1 

260 85 FORMAT ( 4X, D 1 5. 8,X, D 1 3. 6, X, D13 . 6,X, D9 . 2) 

270 WRITE <9,3 5)T9,PI,V9, Dl 

280 T1=T+E1/0.4858D 03 

290 T2=T+E2/0. 1012D 04 

300 T3=T+E3/0.4048D 03 

310 T4=T+E4/0.4858D 03 

320 H=(H+H1* 1D-06)*0.254D 01 

330 P=-C*H*(T2-0.2D 02) + 0. 227D-07*H*H-0. 1 D- 03*H 

340 H=H+P+0.6085D-03*CCT3+T4)/0.2D 1-T 1 ) + • 1 1 08 D- 04* C T 1- • 2D 02) 

350 H=H+0. 1289D-03+0. 1 324503D-05*D 1 

3 55 C6=C1 + C2*T9 + C3*T9**2+C4*T9**3+C5*T9**4 

360 Y=(C0.1D 01+C6*T9)**3)-0. ID 01 

370 W=PI+Y+PI*Y+CT9 + 0.27315D 03+0. 183 053D- 02*H) *V9/V 

38 Z=H+H*W 

390 WRITE C9,240)S0, SliZ 
400 240 FORMATC 'BOOKMX, f STATE' 18X, f 7: f /F6.3,X, F6.2,4X, D17. 10) 

410 GO TO 60 

420 END 

DATA3 

10 10.036 15.01 -0.47D 02 -0.20D 01 -0.15D 02 -0.24D 02 

12 0.199761D 02 0.248000D 02 0. 160D 02 

14 0.32732650D 03 0.629742D-05 0.701470D 00 0. 18D 02 

16 10.036 16.01 -0.50D 02 -0.40D 01 -0.16D 02 -0.25D 02 

18 0.199887D 02 0.209000D 02 0. 176D 02 

20 0.23146721D 03 0.714875D-05 0.786947D 00 -0.50D 01 

22 10.036 16.02 -0.50D 02 -0.60D 01 -0. 15D 02 -0.25D 02 

24 0.199878D 02 0.209000D 02 0. 176D 02 

26 0.23146773D 03 0.714704D-05 0.786771D 00 0.29D 02 

28 10.036 17.01 -0.51D 02 -0.20D 01 -0.16D 02 -0.27D 02 

30 0.199772D 02 0.170000D 02 0.212D 02 

32 0.13616997D 03 0.825908D-05 0.893907D 00 -0.20D 01 

34 10.036 1.04 -0.53D 02 0.00D 00 -0.16D 02 -0.27D 02 

36 0.199949D 02 0.113900D 02 0.960D 01 

36 0.27010000D-01 0.108452D-04 0.112928D 01 -0.10D 01 
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5 DIM ^C90),VC6,90) 

6 DIM EC9 0),U<9 0) 

7 FILES LSDAT 

10 PRINT "INPUT: NO* OF CONSTANTS 1 
20 INPUT M 

110 1=1 

11 1 GO TO 1 15 

112 1=1+1 

115 PEAD #1,BCI ),v(i) 

116 IF END #1 THEN 260 

117 E(I)=E(I)+273. 15 
120 W(I)=1 

130 V( 1*1 )= 1/BCI ) t2 
140 VC2, I)=1/BCI) 
150 V<3*I>*1 
160 U(4,I)=E<I ) 
170 V(5*I )=BCI) t2 
180 V(6*I )=B<I) r3 
250 GO TO 1 12 

260 N=I- 1 

261 PPINT 

262 PFINT "N = M N 

263 PPINT 

27 GO SUE 3 13 

280STOP 

313LETN1 = 

315F0PI=1T0N 

316LETVCM+ 1*I)=YCI> 

319LETN1=N1 + 1 

321NEXTI 

333F0PI=N+ 1TON+M 

344F0PK=1T0M+ 1 

355LETVCK* I )=0 

366NEXTK 

377LETV<I-N*I )= 1 

399NEXTI 

41 1F0PK= 1TOM+ 1 

412IFK=M+ 1THEN417 

415LETT=1 

416G0T0422 

417LET^=- 1 

413G0T0433 

420LETT=-T 

4221 FK= 1 THEN 5 1 1 

433F0FL= 1TOK- 1 

444LETS(L)=0 

455F0PJ= 1TON 

466 S ( L) = S ( L ) ♦ U( Kj J ) * V< L> J ) *V< J ) 
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47 7N EXT J 
488NEXTL 
51 1F0PI= 1T0M+M 
533IFK* 1THEN577 
544F0PL=1T0K-1 

555LETV<K*I)=V<K*I)-S<L>*VCL*I> 
566NEXTL 
577NEXTI 
600IFT<0THEN699 
611LETD(K)=0 
622F0RI=1T0N 

633 D<K)=DCK)+VCK, I >*V(K.> I >*W< I ) 
644NEXTI 

655LETD=SQP(DCK>) 
666F0PI= 1T0N+M 
677LETV(K*I )=V(K^ I)/D 
688NEXTI 
69 3IFT>0THEN420 
699NEXTK 
922LETS=0 
933F0PI=1T0N 

944 S-S+VCM+ 1*1 )*V(M+1^I)*V(I ) 
955NEXTI 

966LETS=SQPCS/(N1-M)) 

1 1 1 1PPINT"TERM", "COEFFI CI ENT"> "STD. DEV. M 
1122PPINT 
1133F0PL= 1TOM 
1144 LET C(L)=0 
1155 FOPK=LTOM 

1166LETCCL)=CCL)+V(K*N+L>*V<K,N+L) 
1177NEXTK 

1188LETCCL)=S*SGP<C(L>> 
1199PRINTL*-VCM+ 1*N+L)*CCL> 
121 1NEXTL 
1222PPINT 
1233PPINT 

1244PPINT"STAMDAPD DEVIATION'S 
1255PPINT 
1266PPINT 

1277 PRINT"B< I ) "* "YC I ) "* "Y f, j "DEV* M 
1288PPIMT 
1299F0RI= 1TON 

1311 ppimt eci>;yci); v ci)-vcm+ ui>; v<m+i,I) 

1313 X=X+V(M+ 1*1) 
1322NEXTI 

1323 PRIN^ 

1324 PPINT "SUM OF RESIDUALS =" X 
13 33 RETURN 

60 00 END 
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1400 457.29, . 079 0,4 57.29/ .08 03, 4 57 . 29, . 08 00, 4 57 . 29, . 08 1 

14 10 457.29, . 076 7,4 57.29, .079 2, 4 57.29, . 0789,4 5 7.29, .0759 

1420 457.29, . 08 16,457. 29, .0796 

14 30 424.98, . 0654,4 24.9 8, . 068 7, 424.9 8, . 063 2 

144 40 0.08, .0607,400.08, . 062 6, 4 00. 08, . 06 1 7, 400. 08, . 62 1 

14 5 39 1.76, .059 6,39 1.7 6, .0 588, 39 1.7 6, .060 1 

1460 379. 56, .0555,379.56, .0567 

147 359. 56, . 0514, 3 59.56, .0507, 3 59. 56, .0528 

148 32 7.32, .049 0, 327. 32, . 049 8, 327 . 32, .049 4, 3 2 7.32, .0479 
1490 323.88, .0431,323-88, .0470 

1500 295.30, .0453, 295.30, .0454 
1510 263. 35, .0457,263.35, .0468 
1520 248. 09, . 0424, 248.09, .0447 
1530 23 1.4 5, . 0443, 23 1 . 4 5, • 046 1, 23 1 . 4 5, . 04 54 

154 199.6 3, .0417, 199.6 3, .04 20, 199.6 3, .04 23 

155 17 2.65, .0380, 172.65, .0403 
1560 167.83, .0389 

157 14 1. 59, . 0372, 14 1.59, .03 54 

158 120.78, .0320, 120.73, .029 6 

159 100, .0247, 100, .0 252, 100, .0 249, 10 0, .0 263, 10 0, .022 0, 100, . 02 3 5 
1600 100, .0274 

1610 3 2.79, .0217,82.79, .0220,82.79, .0178 

1620 62.06, .01 13,62.06, .0 124, 62.06, .0 125 

1630 41.36, .0065,41.36, .0082 

1640 20. 66, .0025, 20.66, .0030 

1650 2.4079, -.00003 

1660 .0886, -.00017 

1670 .2129, .00009 

1680 .0521, .00002 

1690 0.0100,0 
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TEPM 



COEFFICIENT 



STD.DEV. 



1 
2 
3 
4 
5 



120387.7838 
- 1213.5329499 
4.3159552 
-.64407564676E-02 
•35663884587E-05 



669 1.409 68 38 

62.659 08 179 

.2 1 163 09 5 

30556699962E-03 

15950626366E-06 



STANDAFD DEVIATIOM= . 1 52308478 63E- 02 



EC I) 



730.44 
730.44 
730.44 
730.44 
730.44 
730.44 
730.44 
730.44 
730.44 
730.44 
698. 13 



698. 
698. 
673- 
673. 



13 
13 
23 

23 



673.23 
673.23 
664.91 
664.91 
664.9 1 
652.71 
652.71 
632.71 
632.71 
632.71 
600.47 
600.47 
600.47 
600.47 
597.03 



.79000 
.80300 
.80000 
.8 0100 
.76700 
.79200 
.739 00 
.75900 
.8 1600 
.79600 
.65400 
.68700 
.68200 
.6 07 00 
.62600 
.61700 
.62100 
.59600 
.58800 
.60100 
.55500 
.56700 
.51400 
.50700 
.528 00 
.49000 
.498 00 
.49400 
.479 00 
.43100 



Y C I ) 

000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-0 1 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 
000000E-01 



•79392362821E-01 
.79392362821E-01 
.79392362821E-01 
• 79392362821E-01 
•79392362821E-01 
.79392362821E-01 
.79 39 23628 2 IE- 01 
•79392362821E-01 
.79392362821E-01 
.79392362821E-01 
.67446957072E-01 
.67446957072E-01 
•67446957072E-01 
.60436414255E-01 
.60436414255E-01 
.60436414255E-01 
.60436414255E-01 
.58479 128545E-01 
.58479 128545E-01 
.58479 128545E-01 
•55930437313E-01 
.55930437313E-01 
.5251 1784289E-01 
.5251 1784289E-01 
.5251 1784289E-01 
.48688352293E-01 
.48688352293E-01 
.4868S352293E-01 
.48688352293E-01 
.48383358275E-01 



DEV. 

.39236282133E-03 

• 90763717880E-03 
.60763717865E-03 
.70763717862E-03 
•26923628212E-02 
. 19236282138E-03 
•492362S2130E-03 
.34923628212E-02 
. 2207637178 7E-02 
.20763717875E-03 
.20469570721E-02 
. 12530429 280E-02 

• 75304292793E-03 
.26358574537E-03 
.21635357454E-02 
. 12635857453E-02 
. 16635857454E-02 
. 1 1208714551E-02 
.32087145507E-03 
. 16208714551E-02 
.43043731304E-03 
. 76956268 638E-03 
.111 17842889E-02 
. 18 117842888E-02 
.28821571 121E-03 
.31 164770716E-03 
.111 16477072E-02 
.71 164770717E-03 

• 78335229279E-03 
.52833582746E-02 
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597.03 .47000000000E-01 . 48383358275E-0 1 - . 1 3833582745E- 02 
568.45 .45300000000E-01 . 4643 19 548 63E-0 1 - . 1 13 19 548632E- 02 
568.45 .45400000000E-01 . 4643 19 548 63E-0 1 - . 1 03 19 548632E-02 

536.5 .45700000000E-01 . 45060865 144E-0 1 . 639 1 3485577E- 03 

536.5 .46800000000E-01 . 450608 65 144E-0 1 . 1739 1348557E-02 
521.24 .42400000000E-01 . 445 1 1 634947E-0 1 - . 2 1 1 1 6349463E-02 
521.24 .44700000000E-01 . 445 1 1 634947E-0 1 . 183 365053 1 5E- 03 

504.6 .44300000000E-01 . 4386209 679 IE- 1 . 4379 03209 1 6E- 03 
504.6 .46100000000E-01 .4386209679 1E-01 . 22379 03209 2E- 02 
504.6 .45400000000E-01 . 4386209679 1E-01 . 1 537903209 IE- 02 
472.78 .41700000000E-01 . 42088479685E-0 1 - . 38847968548E-03 
472.78 .42000000000E-01 . 4208847968 5E- 1 - .88479 685445E- 04 
472.78 .42300000000E-0} . 42088479685E-0 1 . 2 1 1 5203 1448E-03 
445.8 .38000000000E-01 . 39 57325489 3E-0 1 -. 1 57325489 30E- 02 
445.8 .40300000000E-01 . 39 57325489 3E-0 1 . 726745 1 0700E- 03 
440.98 .38900000000E-01 . 38990228025E-0 1 - .9 0228 024703E-04 
414.74 .37200000000E-01 . 349 5687 19 26E-0 1 . 2243 1 28074 1E-02 
414.74 .35400000000E-01 . 349 5687 19 26E-0 1 .443 128074 1 0E-03 
393.93 .32000000000E-01 . 3 06 1 528 3 1 23E-0 1 . 1 3847 1 68 776E- 02 

393.93 .29600000000E-01 . 306 1 528 3 1 23E-0 1 - . 1 1 5283 1 225E- 02 
373.15 .24700000000E-01 . 252339 18747E-0 1 - . 5339 1874735E-03 
373.15 .25200000000E-01 . 252339 18747E-0 1 -. 339 1874736 IE- 04 
373.15 .24900000000E-01 . 252339 18747E-0 1 - . 3339 18 74734E- 03 
373.15 .26300000000E-01 . 2 52339 1S747E-0 1 . 1 06608 1 2526E-02 
373.15 .22000000000E-01 . 252339 18747E-0 1 -. 32339 187474E-02 
373.15 .23500000000E-01 . 2 52339 18 747E-0 1 - . 1 7339 18 7473E-02 

373.15 .27400000000E-01 . 252339 13747E-0 1 . 2 1 6603 12527E-02 

355.94 .21700000000E-01 . 2007 1 5 1 1472E-0 1 . 1 628483 528 1E-02 
355.94 .22000000000E-01 . 2007 1 5 1 1472E-0 1 . 1928433 5230E-02 
355.94 . 17300000000E-01 . 2007 1 5 1 1472E-0 1 - . 227 1 5 1 147 19E- 02 
335.21 . 1 1300000000E-01 . 1 33 14065430E-0 1 - . 20 140654299E-02 
335.21 . 12400000000E-01 . 1 33 14065430E-0 1 - .9 140654299 0E-03 
335.21 . 12500000000E-01 . 1 33 14065430E-0 1 - .3 1406542992E- 03 
314.51 .65000000000E-02 . 6679 0537379 E-02 - . 1 79 05373787E- 03 
314.51 .32000000000E-02 . 6679 0537378E-02 . 1 5209462622E- 02 
293.81 . 25000000000E-02 . 1 52 1 09896 12E-02 .9 789 1 03884E-03 
293.8 1 .30000000000E-02 . 1 52 1 09896 12E-02 . 1 4739 1 0338E- 02 
275.5579 - . 30000000000E- 04 .9 46 1 6373689E- 04 -. 1 246 1 637369 E- 03 
273.2336 - . 1 7000000000E-03 . 25398 657250E- 03 - .42398 657250E- 03 
273.3629 .9 0000000000E-04 . 2430 1 654224E- 03 - . 1 530 1 654224E-03 
273.2021 .20000000000E-04 . 25726 1 503 10E- 03 -. 23726 1 503 1 0E- 03 

273.16 .261069 13775E-03 -. 26 1 069 1 3775E-03 

SUM OF PESI DUALS = -. 4774847 1843E- 1 1 
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The least squares solution presented does not give a calculated 
value of zero for T 68 /K 68 - T/K at T 68 = 273. 16 K 6Q , but amounts to 
0*000261. The equation has a minimum at 278*57 K/-o* where the value 
is less than 0*000023. However, the variation at 273»l6 K^n is 
within experimental error, and can be ignored for any scale adjustments. 
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